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Abstract We consider smoothing equations of the form 

x ^ j2 T i x i+ c 

3 > 1 

where ( C , Tf, X 2 ,...) is a given sequence of random variables and Xi . X2 ,... 
are independent copies of X and independent of the sequence (C,Ti,T %,...). 
The focus is on complex smoothing equations, i.e., the case where the random 
variables X, C, Ti, T 2 ,... are complex-valued, but also more general multivari¬ 
ate smoothing equations are considered, in which the Tj are similarity matrices. 
Under mild assumptions on (C,Ti,T 2 ,...), we describe the laws of all random 
variables X solving the above smoothing equation. These are the distributions 
of randomly shifted and stopped Levy processes satisfying a certain invariance 
property called (U, a)-stability, which is related to operator (semi)stability. 
The results are applied to various examples from applied probability and sta¬ 
tistical physics. 
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1 Introduction and main results 

1.1 Smoothing equations 

In the paper at hand, we consider complex smoothing equations of the form 

x ^J2 t j x j+ c ( 1 - 1 ) 

3> 1 

where '=" denotes equality in law, (C. T\. T ?,...) is a given sequence of com¬ 
plex random variables and X\,X 2 , ... are independent copies of the complex 
random variable X and independent of ((7, 7"i, T 2 ,...). The law of X is called 
a solution 1 to (1.1). Under suitable assumptions on (C, Tf, X 2 ,...), we provide 
a complete description of the set of all solutions to (1.1). 

This kind of problem has a long history going back at least to [35,41,50]. 
Recently, there has been progress leading to a complete description of the set 
of all solutions to (1.1) in the case where all the Tj, j £ N are real and C and 
X are random vectors [5,43]. We also refer to [2] for an overview over the early 
results in the field. 

The complex case can naturally be embedded into the more general mul¬ 
tivariate case, in which C and X are d-dimensional and the Tj are random 
d x d matrices, d £ N. In this setup, (1.1) has been solved under two different 
sets of assumptions, namely, in the homogeneous case (i.e., (7 = 0 a.s.) under 
assumptions that guarantee that all solutions are essentially scale mixtures of 
multivariate normal distributions [14], and in the case that (7,71,72,... have 
positive entries only and attention is restricted to solutions X on [0, oo) d [64]. 
Yet, these results either cannot be applied to the complex case [64] or cover 
only a very special situation [14], We will solve (1.1) in a multivariate setup 
that comfortably covers the complex case and thereby address open questions 
in papers by Barral [10, Remark 4] (posed for (7 = 0 and 7\, T 2 ,... taking 
values in a Banach algebra); Chauvin et al. [29, Remark 4.5] and Madaule et 
al. [61, Section 1.2] (in the complex case). 


1.2 Motivation 

We believe that solving (1.1) in a setup as general as possible is of high the¬ 
oretical value for probability theory. Smoothing equations appear naturally 
in various fields of probability and statistical physics. We mention here two 
classes of examples that will also motivate the choice of our setup. 

1.2.1 Complex smoothing equations in models of Applied Probability 

In various models of Applied Probability including 6-ary search trees [29,30,37, 
46], Polya urn models [46,53,69], the conservative fragmentation model [47], 

1 In slight abuse of language, we will sometimes call a random variable X a solution if 
the distribution of X is a solution to (1.1). 
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and B-urns [28] phase transitions were shown for the limiting behavior of the 
distributions of quantities of interest. All these phase transitions follow very 
similar patterns suggesting that they are particular instances of one universal 
phenomenon. 

Before we sketch the general phenomenon, we describe it in the context of 
b- ary search trees. The space requirement (i.e., the total number of nodes) in 
a b- ary search tree with n keys inserted under the random permutation model, 
centered by its mean and scaled by y/n, is known [58] to be asymptotically 
normal when b < 26. On the other hand, it exhibits stable periodic fluctuations 
[30,37] around its mean when b > 26. 

The general scheme is as follows. For each of the models mentioned above, 
a characteristic equation can be formulated with several complex roots (for 
6-ary search trees, b appears as a parameter in the equation). The root with 
largest real part always is 1. The asymptotic behavior of the quantity of interest 
is determined by the root with second largest real part, £ + i rj, say. Roughly 
speaking, if £ < i, the fluctuations around the expected size of the quantity are 
of order y/n. (times possibly a slowly varying correction term), where n denotes 
the “size” of the model in an appropriate sense, and the limiting distribution 
is normal. When £ > i, the fluctuations around the mean are of order rfi and 
there is no convergence but a periodic limiting behavior the precise description 
of which involves solutions to complex smoothing equations. Our results reflect 
the phase transition between normal and periodic limiting behavior with stable 
fluctuations. Concrete examples are considered in detail in Section 2. 

1.2.2 Kinetic models 

A main objective in the kinetic theory of gas is to understand the distribution 
[it on R 3 of particle velocities and its evolution in time, given by the equation 

—[it + [it = S([i t ). (1.2) 

Eq. (1.2) has to be understood in the weak sense, i.e., 

^ J fd(H + J fd/J-t = j fdS(nt) 

for all bounded continuous functions / on R 3 . Here S' is a mapping on proba¬ 
bility measures, called collisional-gain operator, that describes the change in nt 
due to the collision of two uncorrelated particles. Its form can be deduced from 
the equations of energy and momentum conservation (see [73] for a detailed 
survey on the topic). 

In order to understand qualitative aspects of (1.2), simplifications of S 
have been considered, the most important one being the Kac caricature of the 
Boltzmann equation [48], where [it is restricted to the real line, and 

S : [j. i—^ law of (sin(<9)Vj + cos(<9)1-2) 
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where Vj, Vi are i.i.d. with law p and assumed to be independent of the random 
arc length O which has the uniform distribution on [0,27 t). In this case, only 
energy is conserved, but not momentum. Any steady state distribution Hoc, 
i.e., satisfying = 0, is a solution of the smoothing equation S(p oo) = Moo 

with real valued weights. Such equations were investigated in [43]. 

Recently, there has been a lot of interest in generalizations of the Kac 
caricature [11,12,22], also with applications in economic theory [27,32,63], to 
mention just a few. 

Of particular interest for us is a generalization of the Kac model in R 3 , 
considered by Bassetti and Matthes [14, Section 6.2]. Here, a random vector 
V the law of which is the steady state distribution, satisfies 

V '= LV 1 +RV 2 , (1.3) 

where V. V\, Vi are i.i.d. and independent of the random pair (L, R ) of similar¬ 
ity matrices 2 which is subject to the condition E[||L|| +||i?||“] = 1 and density 
assumptions, which guarantee that there is a unique solution (up to scaling) to 
(1.3). It is shown that this solution is a mixture of multivariate Gaussian dis¬ 
tributions. We solve (1.3) in much greater generality. In particular, dropping 
their density assumption, we show that larger classes of Gaussian solutions 
appear. 


1.3 Setup and main results 

We aim for solving (1.1) in a multivariate setup covering both, the complex 
equation and the multivariate equation as outlined in Section 1.2.2. 

Fix d £ N, denote by |x| = ( x\ + ...a^) 1 / 2 for x = (x\,.. ., Xd) £ R d 
the Euclidean norm and by ||a|| := sup{|a:r| : |a:| = 1} for a real d x rf- 
matrix a the associated matrix norm. Assume that we are given a sequence 
(C,Ti,Ti ,...) =: ( C,T ) where C is a d-dimensional random vector and T = 
(T)) je N is a sequence of random d x d similarity matrices, i.e., for each j £ N, 
Tj = ||Tjj| Oj where Oj £ O(d), with O (d) denoting the group of orthogonal 
dx d matrices. Notice that ||T)|| and Oj are random elements of R> = [0, oo) 
and O(d), respectively, which in general depend on each other. Further, we 
point out that no assumptions on the dependence structure of the sequence 
(C, T ) are imposed. 

This setup includes the case where X, C, Xj, Ti,... are complex random 
variables, henceforth referred to as the complex case, since C can be identified 
with R 2 and multiplication by a complex number re lS , r > 0, 9 £ [0, 27r) 
corresponds to multiplication in R 2 from the left by the similarity matrix 

/ rcosO —r sind\ 
yrsind rcosO J 

2 By a d X d similarity matrix we mean a d X d, matrix that can be written as a scale 
multiple of an orthogonal matrix. 
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We proceed by introducing the assumptions and some basic concepts needed 
for the statement of our main result. 


1.3.1 Assumptions 


Throughout the paper, we assume that the number of non-zero Tj is a.s. 
finite, i.e., N := #{j € N : ||Tj|| > 0} < oo a.s. We suppose without loss of 
generality that Ti,..., T/v ^ 0 and Tj = 0 for j > N. For the formulation of 
further assumptions, we define, for s > 0, 


m(s) 


E 


N 

Xjia* 

3 =1 


Tj ir 


(1.4) 


Notice that m(0) = E[7V]. Throughout the paper the following assumptions 
will be in force. 


E [N] > 1. (Al) 

m(a) = 1 for some a > 0. (A2) 

(Notice that we do not exclude the case that m(s) = oo for all s ^ a.) Then 
W\ := Y!*_ x \\T.j ||° is a nonnegative random variable with unit mean. In our 
main results, we assume that 


m {a) 


:=E 


N 

Elixir ^ n T iii 


i=i 


G (—oo,0) and E[VFi log + W\) < oo. 


(A3) 


Notice that m'{a) as defined in (A3) is indeed the derivative of s i-)- m{s) at 
a if the latter exists, a can be thought of as a generalized index of stability. 
Indeed, strictly a-stable random variables solve the particular instance of (1.1) 
with C = 0 and Tj = N -1 / a , N > 2 fixed. 





(a) m(s) for Eqs. (2.3) and (b) m(s) for Eq. (2.8) with 
(2.4) with b = 27 b = 7 


(c) m(s) may eventually in¬ 
crease if sup j \Tj\ > 1 with 
positive probability 


Fig. 1.1: Possible shapes of s i-> m(s) 
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While the assumptions imposed up to now are very mild and suffice when 
a ^ 1 , it might be not surprising that additional assumptions are needed in 
the case a = 1 in order to overcome severe technical obstacles. We recommend 
to skip the following part at first reading. 


1.3.2 Additional assumptions for the case a = 1 

We solve (1.1) in the case a = 1 under two sets of assumptions, which cover 
the most relevant cases. Let £ = 1 + diniR0(d) = 1 + d(d — l)/2. Then the 
first set of additional assumptions for a = 1 is 


E 


N 




L i=i 


is spread-out (w.r.t. the Haar measure on S(d)), 


r N 


E 


Eii^iriio g -(ii7) 


\t+s+l 


3 =1 


< oo and E[/i 2 £f$+i(Wi)] < oo 


(A4) 


for some S > 0 and h r (x) := ^log 4 ’^))’’log + (log + (:r)). Notice that in the 
complex case, for the first condition in (A4) to hold, it is sufficient that, for 
some j £ N, the law of Tj is spread-out on C. Let O denote the smallest closed 
subgroup of O (d) which contains the random set {Tj/ ||Tj|| : j = 1,..., N} 
with probability one. Assumption (A4) will imply that O 3 SO (d), the sub¬ 
group ofO(d) of matrices with determinant 1. In the second set of assumptions, 
we assume O to be a finite group: 


is finite, E 


N 


E ll T fir <5 -i°g(ll'L-|l)( 

L i=i 


r N 


E 


Elixir (log-diT.i 

3 =1 


is spread-out, 

..... j 

21 < oo and E[/i 3 (Wi)] < oo. 


(A4’) 


For a probability measure /i on O (d) x R, we say that p. satisfies the minoriza- 
tion condition (M) if there is a nonempty open interval / C R and 7 > 0 such 
that 

p(do,dx) > 7lso(d)x/(°i x ) H@( d )(do) da; (M) 

where ILo(d) is tl ie normalized Haar measure on the compact group O(d). We 
will show in Lemma 3.6 that we may assume the following stronger property 
instead of (A4): 


r N 

E ||,— log||Tj II) (*) 

l j =1 


(A4) holds and E 


satisfies (M). (A5) 
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1.3.3 Weighted branching 

The weighted branching process is a natural tool in the study of equations of 
the form (1.1). In our context, this process is defined as follows. 

Let V = UneNo denote the infinite Ulam-Harris tree where N° := {0} 
is the set that contains the empty tuple only. For u,v G V, u = (iq,..., u m ), 
v = (vi,..., v n ), we write uv for (u±, ..., u m , Vi,..., v n ). We say that v is in 
generation n, in short: |u| = n, if v G N™. The restriction of v to its first k 
components is denoted by v\k- 

Let (C, T) := ((C{v),T(v))) v&f = ((C(v),Ti (v),T 2 (v), .. .))„ev be a fam¬ 
ily of i.i.d. copies of the sequence (C,T). For notational simplicity, we as¬ 
sume (C(0),T(0)) = ( C,T ). Let Id denote the dx d identity matrix and let 
L(0) = Id- Recursively, for v G V and j G N, we define L(vj) = L(v)Tj(v). 
Hence, if v = v\ ... v n G N", then 

L(v) = T Vl (0) ■ ... ■ T Vn (v | n -i). 

Notice that the order of multiplication matters. Each of the matrices Tj(y) 
and L(v) (if nonzero) can be written as the product of a positive scaling factor 
and an orthogonal matrix: 

Tj(v) = e~ Sj ^Oj(v) and L(v) = \\L(v)\\0(v) 

where 

Sj(v) = - log ||T j (u)|| € R, Oj(v) = ||T j (u)||- 1 T i (u) G O(d) 

S(v) = - log ||L(u)|| G K, 0(v) = HL^ir 1 L(v) G O (d) 

whenever ||Tj(u)|| > 0 or ||T(u)|| > 0, respectively. We make the convention 
that whenever we quantify over the |u| = n as in or ll| v |=n’ ^ ias 

to be understood as a quantification over the |i>| =n with ||L(u)|| > 0 only. 
These definitions imply that, for v = iq ... v n , when ||L(u)|| > 0, 

n 

S{v) = y ^S Vk (i^lfc-i) and 0(v) = O Vl (0) •... • 0 Vn (v | n -i)- 

k =1 

Here, of course, the order of summation does not matter while the order of 
matrix multiplication does (in general). We further point out that conditions 
(A1)-(A2) imply that 


lim 

n—^oo 


sup ||L(u)|| 

|h| =71 


0 a.s. 


(1.5) 


(with the convention that sup0 = 0), see [19, Theorem 3] for a reference. 

For u G Y and a function = ^((C, T)) of the weighted branching process, 
let \&] u be defined as 'P(((C{uv),T(uv))) V £y), that is, the same function but 
applied to the weighted branching process rooted at u. The [•]„, u G Y are 
called shift operators. 
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1.3.4 Special solutions to smoothing equations 

The solutions to (1.1) are connected with the solutions of two related distri¬ 
butional identities, namely, the tilted homogeneous equation for nonnegative 
random variables 

w = w Elixir w i (1-6) 

3 > 1 

where W\ 1 W 2 , ■ ■ . are i.i.d. copies of the nonnegative random variable W and 
independent of T, and the homogeneous equation 

X ^ ^ 3 X 3 (1-7) 

3 > 1 

where X\,X 2 , ■ ■ . are i.i.d. copies of X and independent of T. Special solutions 
to (1.6) and (1.7) can be constructed using the weighted branching process. 
Indeed, (A2) implies that 

W n := £ ||L(i>)|| a = ]T n € N 0 (1.8) 

\v\=n \v\=n 

defines a nonnegative mean-one martingale. Let W = limn^oo W n a.s. It is 
known [60] that (A1)-(A3) guarantee E [W] = 1. It can be checked that 

W= \\mt[W} v a.s. (1.9) 

\v\=n 


for every n G No- In particular, IT is a solution to (1.6). The set of solutions to 
(1.6) is (Law(cIT) : c > 0}, see [2,42]. The description of the set of solutions 
to (1.7) is more delicate and most of the analysis in this paper is concerned 
with solving it. One aspect of this equation is that special solutions may arise 
due to balancing effects in the sum Z\ := >1 Tj. This may happen when the 

matrix E[i?i] has eigenvalue 1. Then, for any eigenvector w corresponding to 
the eigenvalue 1 and with Z n := L(v), the sequence (Z n w) n ^ 0 defines 

a martingale. If it converges in probability, we denote its limit by Z w . Z w is a 
function of T and a solution to the identity 

Z = a.s. (1.10) 

3 > 1 

The solutions to (1.10) are described by the following result. Note that Z = 0 
a.s. always satisfies ( 1 . 10 ), this we call the trivial solution. 

Proposition 1.1 Assume that (A1)-(A3) hold and let Z satisfy (1.10). 

(a) If 0 < a < 1, then Z = 0 a.s. 

(b) If a = 1 and (A4) holds, then Z = 0 a.s. If (A4’) holds, then a nontriv¬ 
ial solution Z exists iff w — E [Z] is an eigenvector corresponding to the 
eigenvalue 1 ofE[Zi\, and then Z = Ww. 
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(c) If l<a<2, then a nontrivial solution Z exists iff w = E [Z] exists and is 
an eigenvector corresponding to the eigenvalue 1 o/E[Zi ] and (Z„tt>) ra€ N 0 
is uniformly integrable. If these conditions hold, then Z = Z w and the 
sequence (ZnU^jigNg is bounded in C s for all 1 < s < a. 

If 1 is an eigenvalue of E[Zi] associated with the eigenvector w, then a 
sufficient condition for CJ 5 -boundedness (and hence uniform integrability) 
of (Z n w) n £Fj 0 for ft € (a, 2 ] is E[|Ziu>|^] < oo and m(/3) < 1. 

(d) If a > 2, then a nontrivial solution Z exists iff there is a deterministic 
w ^ 0 such that w is an eigenvector corresponding to the eigenvalue 1 of 
the random matrix Z\ a.s., and then Z = w a.s. 

Besides W and Z, the class of (U, a)-stable Levy processes, with U being a 
closed subgroup of the group of similarity matrices of R d , will be relevant for 
the description of solutions, and is introduced next. This notion is a particular 
case of ([/, (instability, a concept introduced in [23, p.338]. 

Let (l*)t>o be a Levy process on R d . We say that (Y t )t> o is ( U, a)-stable 
if there exists a mapping b : U —> R d such that, for all u £ U, t > 0, 

uY t 'Iff Y lHrt +tb(u). (1.11) 

This implies in particular that each Y t is operator semistable [39], see Section 
4.2 for more details. Property (1.11) is equivalent to 

I / {u T x) = ||w|| Q \h(x) + i{x, b(u)) (1.12) 

for all u € U for the characteristic exponent W of Yi. (Yt ) t > 0 is said to be 
strictly (U, a)-stable if b(u) = 0. We say that a probability measure P on R d is 
(strictly) (U, a)-stable if it is the law of a (strictly) (U, a)-stable Levy process 
at time 1 . 

The general formula of the Levy measure of a (U, a)-stable law is given in 
Proposition 1.6 below. Let us just point out that there are no (U, a)-stable 
Levy processes for a > 2 unless U C O(d) and that the larger the group 
U, the smaller the class of (U, a)-stable Levy processes. For example, if U D 
R> x SO (d), then \P(x) = —c||x|| Q for some c > 0, i.e. (Y t ) t >o is strictly a- 
stable and rotation-invariant. It will be shown (in Remark 3.7 below) that 
assumption (A4) implies U D R> x SO (d). Below, we consider (U, a)-stable 
Levy processes, where U denotes the smallest closed subgroup of the similarity 
group that covers {Tj : j = 1,..., N} with probability one. 

There are further technicalities to deal with before the main result can 
be formulated but the complex and homogeneous case, the most important 
special case in view of applications, is now given as an illustration. 

1.3.5 Solutions to complex smoothing equations 

In the complex case, each Z n is a complex random variable, and the condi¬ 
tion that E[Zi] has eigenvalue 1 is equivalent to E[Zi] = 1, which in turn is 
equivalent to (Z n ) ra€ N 0 being a (complex) martingale. We write Z for the a.s. 
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limit of this martingale if it exists, and define Z = 0, otherwise. Note that U 
is now the smallest closed subgroup of the multiplicative group C* that covers 
{Tj : j = 1,... ,N} with probability one. 

Theorem 1.2 Consider (1.7) in the complex case and assume that (Al)-(A3) 
hold. Additionally suppose that (A4) or (A4’) holds if a = 1. 

Then a probability distribution on C is a solution to (1.7) if and only if it 
is the law of a random variable of the form 

Y w + aZ (1.13) 

for some a £ C and a complex strictly (U, a)-stable Levy process (Y)t>o inde¬ 
pendent of (W , Z). 

Theorem 1.2 is a special case of the more general Theorem 1.5. Therefore, 
we do not give a separate proof of Theorem 1.2. 

We have tried to present the result of Theorem 1.2 as concise as possible. 
This should not hide the fact that it hosts a multitude of different cases. On 
the one hand, there are several qualitatively different possibilities for the group 
U two of which, coming from examples discussed in Section 2.1, 



(a) The group U associated with Eq. (2.4) (b) The group U associated with Eq. (2.8) 

appearing in the context of b -ary search appearing in the context of cyclic Polya 
trees (here, b = 27) urns (with b = 7) 

Fig. 1.2: The group U in examples 

The reader should notice that U may also be a family of discrete points aligned 
on a “snail graph” or consist of a finite number of copies of R> := ( 0 ,oo) 
obtained by multiplication with roots of unity as depicted in Figure 1.3. On 
the other hand, different values of a give rise to qualitatively different regimes. 
A brief discussion of the implications of Proposition 1.1 and the structure of 
(U, a)-stable Levy processes (described in Proposition 1.6 below) is given in 
the next remark. 

Remark 1.3 Before we discuss the different cases of Theorem 1.2, we point out 
that (A2) implies that U % S. 
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(a) Example where HJ is discrete 



(b) Example where U consists of circu¬ 
larly arranged copies of R> 


Fig. 1.3: Two possible shapes for U 


(i) a<l: Then Z = 0 a.s., i.e., the second summand in (1.13) vanishes. 

(ii) a= 1: Z vanishes unless E[Zi] = 1. E[Zi] = 1 and ro( 1) = 1 imply Tj £ 
[0,oo) a.s. for all j £ N. In this case, Z = W and U C R>. In fact, 
Assumptions (A4) and (A4’), resp., imply that U = R>. 

(iii) l<a<2: Z vanishes unless E[Zi] = 1. In most applications, E[Zi] = 1 will 
hold and the full spectrum of solutions given by (1.13) will arise. 

(iv) a = 2: In this case, (Y t )t >o is a centered Gaussian process (or 0). Moreover, 
(U, 2)-stability necessitates that real and imaginary part of (Y t ) t >o are i.i.d. 
centered one-dinrensional Brownian motions (or 0) whenever U \ R ^ 0. 
Further, Z = 1 a.s. iff Z\ = 1 a.s. and Z = 0, otherwise. 

(v) a > 2: Here, Y t = 0 a.s. for all t and Z = 1 a.s. iff Z\ — 1 a.s. and Z = 0, 
otherwise. Hence, the set of solutions to (1.7) is either {d 0 : a £ C} or {<5o}, 
respectively, where here and throughout the paper, S a denotes the Dirac 
distribution with a point at a. 

Remark 1.4 In many cases, the different solutions to (1.13) can be distin¬ 
guished via their tail behavior. We discuss here in details the most relevant 
case 1 < a < 2. In applications, typically E[|Z|^] < oo for some (3 > a (see the 
sufficient condition of Proposition 1.1(c)), whereas Y w for a non-trivial (U, ca¬ 
stable Levy process (Y t ) t > 0 exhibits heavier tails. Indeed, the Levy measure 
v of (Y t ) t >o satisfies h\r~ a < S'({|a;| > r}) < h 2 r~ a for 0 < hi < /12 < 00 , 
which can be derived directly from the (U, a)-stability and is also implicit 
in our proofs. According to [72, Corollary 25.8], this implies that, for every 
t > 0, Yt has all absolute moments of order < a finite, while E[|Yj|“] = 00 . 
Now if G := {|M| : u £ U} = R>, then, using that E[W] = 1 < 00 and the 
independence of W and (l*)t>o, 

E[|ivn = E[(IT 1 /a |Yi|) p ] = E[W p/ “]E[|r 1 | 1 ’] < 00 

iff p < a. In the case where G = r z for some r > 1, one can argue similarly 
using the fact that E[sup 1<s<r |Y S | P ] < 00 iff E[|1^] < 00 , see [72, Theorem 
25.18]. 
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More results on the tail behavior of Z (under stronger conditions than 
imposed here) can be found in [25]. 


1.3.6 Solutions to multivariate smoothing equations 

We continue with the description of the solutions to (1.1) in the general situa¬ 
tion. A special solution can be constructed in terms of the weighted branching 
process. Define 

W: := ^ L(v)C(v), n e N 0 (1.14) 

\v\<n 

and let W* denote the limit in probability as n —> oo of W* provided the limit 
exists. If it exists, W* satisfies 

W* = J2 t AW*}j+ c a.s. (1.15) 

j>i 


and thus constitutes a solution to (1.1). Each of the following is a sufficient 
condition for the convergence in probability of W* taken from [43, Proposition 
2 . 1 ], which remains valid in the present context: 

(51) (Al) and (A2) hold and there is /3 € (0,1] with m(/3) < 1 and IE[|C'|' S ] < oo. 

(52) For some (3 > 1, sup ragN() E^IE*!^] < oo and either Tj > 0 for all j £ N or 
E[C7] = 0. 

Notice that if C = 0 a.s., then W* converges trivially to W* = 0 a.s. 

We now state the main result for the general case. 

Theorem 1.5 Assume that (A1)-(A3) hold and that W* —> W* in probability 
as n —► oo. If a = 1, assume that (A4) or (A4’) holds in addition. 

Then a probability distribution on is a solution to (1.1) if and only if it 
is the law of a random variable of the form 

W*+Y w + Z (1.16) 

where Z is a solution of (1.10) and {Yt)t> o a strictly (HJ ,a)-stable Levy 

process on independent of (W*, W, Z). 

Denote by Ei C the eigenspace corresponding to the eigenvalue 1 of 
the matrix K[Z\] and let E\ = {0} if 1 is not an eigenvalue of E[Zi]. Then 
Proposition 1.1 yields that the set of solutions to (1.10) is either empty or 
parametrized by E\. Besides, the solutions are parametrized by the different 
strictly (U, a)-stable distributions, which always include Y t = 0, i.e., W* + Z 
is also a fixed point. In previous works, the solutions to smoothing equations 
on the nonnegative halfline [4], on R [5] and on in the case where OCR* 
[43] have been represented in the form 


W* +W 1/a Y + Z 


( 1 . 17 ) 
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with (W*,W,Z) defined as here and Y denoting an independent strictly ca¬ 
stable random variable (where Y = 0 is allowed). The same representation 
is possible here when K> x {Id} C U since then (Y t )t>o is strictly a-stable 
(or zero) and thus Y\y has the same law as W x l a Y\. In general, solutions to 
(1.1) do not possess a representation of the form (1.17) as U need not contain 
R> x {Id}, see e.g. the examples depicted in Figure 1.2. 

Observe that W* and Z are measurable functions of (C, T) and satisfy 
a.s. versions of the inhomogeneous or homogeneous fixed point equation, i.e., 
(1.15) and (1.10), respectively. Such fixed points are called endogenous , a no¬ 
tion coined by Aldous and Bandyopadliyay [1], see [2, Section 6] and [43, 
Section 3.5] for further information. The fixed point Y\y is not endogenous, 
for the process (Tt)t>o introduces additional randomness. Note, however, that 
W is an endogenous fixed point of the one-dimensional smoothing transform 
with scalar weights ||Ti||“ ,..., ||Tjv||“, see (1.9). The tail behavior of W* is 
investigated in [25], and parallels that of Z. 


1.4 The class of strictly ( U , a)-stable Levy processes 

To complement our main result, we determine the form of the characteristic 
exponent of strictly ( U , a)-stable Levy processes in this section. The analysis 
of (17, a)-stable Levy processes for U C O (d) is of no relevance for this paper 
and simpler than in the case U ^ O(d) and therefore omitted here. 

Below, denote by i/ a , 0 < a < 2, a Levy measure satisfying 

v a {uB) = \\u\\~ a v a (B) (1.18) 

for all u £ U and all Borel sets B C R d \ {0}. We call such a Levy measure 
(U, a)-invariant as well. The structure of such measures is described in Section 
4.3. Given v a . define the functions 

Vi 0) : = j^j ^ j (i-cos ((x,y)))v a (dy), (1.19) 

V%( x ) '■= J (sin((x,j/)) - l {a>1} (x,y))is a {dy). ( 1 . 20 ) 

It is proved in Lemma B.l that r/,“ are bounded functions, satisfying iji(u T x) = 
r]i(x) for all u € U, x £ \ {0}, i = 1,2. 

Now consider the case where the image of U under the homomorphism 
U 3 « i-> ||u|| is R>. Notice that the groups depicted in Figure 1.2 are of 

this type. It is proved in [24, Proposition C.l and Theorem D.13], see also 

Proposition 4.1, that then U = { t® : t £ R>} x C for a suitable d x d -matrix 
Q and C '■= U (~l O (d). Here, t® := and we scale Q in such a way that 

11| = t. Observe that (1.11) then implies that any ([/, a)-stable law is also 
operator stable with exponent —Q, see Section 4.2 for details. 
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Further, let p be a measure on S d_1 := {x £ R d : |ar| = 1} satisfying 
p{oB) = p(B) for all o £ C and Borel sets B C § d_1 . We call such a measure 
C-invariant and define 


{Q - hW : = LLj fiy) iiw' (d,|ii ' 


( 1 . 21 ) 

( 1 . 22 ) 


given the right-hand side in (1.22) is in the range of (Q — Id). 

Below, we write O := {u/ ||it|| : u £ U} for the projection of U onto O (d). 


Proposition 1.6 Let U %. ©(d) be a closed subgroup of the similarity group 
and a > 0. Then for a d-dimensional Levy process (Ft)t>o with characteristic 
exponent T, the following assertions hold: 

(i) Let 1 a £ (0,2). (Y t ) t >o is strictly (U, a)-stable iff If is of the form 


L>{x) = -\x\ a rj ( f(x) + i\x\ a ^{x) (1.23) 


for a (U, a)-invariant Levy measure v a . 

(ii) Let U = {t® : t e M>} x (7. (Y t ) t > o is strictly (U, 1 )-stable iff T is of the 
form 

L , (x) = 77 1 (ar) + i(y 1 + z, x) (1-24) 

for a C-invariant measure p on S d_1 , satisfying f (x, s)p(ds) = 0 for all x 
with Q T x = x (this also guarantees the existence of j 1 ), and any vector z 
satisfying uz = ||it|| z for all u £ U. 

(iii) (Yt)t >0 is strictly ( U,2)-stable iff'T is of the form 

T(x) = -x T Ex/2, rrer* (1.25) 

for a positive semi-definite symmetric dx d matrix E satisfying oEo T = E 
for all o £ O. 

(iv) Let a > 2 . (lt)t >0 is strictly (U, a)-stable iff \T(x) = 0 for all x £ 
equivalently, Y t = 0 a.s. for all t > 0. 

Notice that for a = 1, the proposition excludes the case where the image of 
U under the homomorphism U 9«i-> ||u|| is a discrete subgroup of R>. This 
is because (A4) and (A4’), assumed in the case a = 1, imply that G = R> 
and hence the discrete case is of no relevance here. 

A description of the matrices E satisfying oEo T = E for all o £ O is 
given in Proposition 4.4. In particular, if there is no proper subspace V C R d , 
satisfying oV = V for all o £ O, then A is a scalar multiple of the identity 
matrix. 
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1.5 Further organization of the paper 

In Section 2, we apply our main results to the examples mentioned in the 
introduction as well as to further important models. The rest of the paper is 
then devoted to the proofs of the main results, with the major, probabilistic 
part of the proof being given in Section 3, at the beginning of which we will 
also introduce further notation and concepts relevant for the proofs. More 
algebraic considerations, concerned with the structure of (U, a)-stable Levy 
processes, are contained in Section 4. The appendix contains a Choquet-Deny 
lemma for functions on U and a rate-of-convergence result for Markov renewal 
processes; the latter result will be needed only in the case a = 1. 


2 Applications of the main results 

In this section, we discuss the examples from the introduction, as well as 
further applications of our main results to the study of Biggins’ martingale 
with complex parameter or Gaussian multiplicative chaos. 


2.1 Applications of Theorem 1.2 

2.1.1 b-ary search trees 

b- ary search trees are 6-ary trees which are basic data structures in computer 
science used in searching and sorting, see [62] for the definition and back¬ 
ground information. Each node of a b- ary search tree can store up to b — 1 
elements from a set of distinct real numbers x ±,..., x n , called the set of keys. 
For the problems considered, it constitutes no loss of generality to assume 
{xi,... ,x n } = {l,...,n}. Denote by Y n the space requirement of a b- ary 
search tree under the random permutation model, i.e., (aq,... ,x n ) is a uni¬ 
form permutation of the set {1,..., n} and Y n is the number of nodes in the 
resulting b- ary search tree. 

Let b > 4. The asymptotic behavior of Y n is coded in the equation 

6-1 

X{z) : = + = °> zeC - ( 2 - 1 ) 

i =i 

The root with largest absolute value is Ai = 1. Let A 2 denote the root with 
second-largest real part and Im(A 2 ) > 0. Then Re(A 2 ) £ (0,1). If Re(A 2 ) < 
1/2, equivalently, b < 26, then after centering and norming, Y n is asymptoti¬ 
cally normal, see [58]. By using martingale methods, Cliauvin and Pouyanne 
[30] have shown that if Re(A 2 ) > 1/2, equivalently, b > 27, then 

Y n = const • n + 2Re(n A2 A) + o(n Re ^ A2 ^) (2.2) 

where o(n Re ^ A2 - ) ) is a term that, after dividing by n Re ^ 2 \ tends to 0 a.s. and 
in £ 2 , and A is a complex-valued random variable. Since n Aa is complex, 
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n Re f A2 )(F„ — const • n) does not converge in distribution, but behaves like 
2Re(?r Im ( A2 )X). 

Fill and Kapur [37] showed that X solves the smoothing equation 

X ‘= Y^Vj^Xj (2.3) 

3=1 

where V\,... ,Vb are the spacings of b — 1 i.i.d. random variables uniformly 
distributed over (0,1), U \,..., £4-i, say. In other words, V) = /7y) — U(j- 1 ) 
for j = 1 ,... ,b with t/( 0 ) = 0, = 1 and ..., U^-i)) being the order 

statistics of (Ui ,..., E4_i). What is more, Fill and Kapur showed that (the 
law of) X is the unique solution to (2.3) subject to the additional constraints 
E[X] = n and E[| AT | 2 ] < oo where 0 is a given complex constant. 

Later, using an embedding into continuous-time multi-type Markov branch¬ 
ing processes, Clrauvin et al. [29] established a connection between the law of 
X and the complex smoothing equation 

X '= e~ X2T (X 1 + ... + X b ) (2.4) 


where T has the same distribution as the sum n + ... + T(,_i of independent 
random variables with tj being exponentially distributed with parameter j, 
j = 1,..., b— 1. Again, the connection concerns the solution to (2.4) with fixed 
expectation E[X] = / 0 and finite second moment E[|X| 2 ] < oo. 

Here we will consider Eq. (2.3) only; the study of (2.4) bears a striking sim¬ 
ilarity. One can check that V[,... .Vi, are identically distributed with Lebesgue 
density (b — 1)(1 — x) b_2 l( 0j i)(a;). Hence, for z £ C, 

E[V? + • • • + ni = Kb - 1) jf 1 *'(1 - x) b ~ 2 dx = bl r { [ Z +Q (2-5) 

6 ! 

(z + 1)- ...-(z + b-1) 

where r denotes Euler’s gamma function. We conclude that in the present 
context the function s <—> m(s) defined in (1.4) takes the form 


m(s) 


HI 


3 = 1 


1 

Re(A 2 )s + j ’ 


( 2 . 6 ) 


Figure 1.1a shows the graph of s <—> m(s). Notice that m(a) = 1 only for 
a = l/Re(A 2 ). In particular, the phase transition at b = 26 is visible here 
since a > 2 for b < 26 and a £ (1, 2) for b > 27, and the transition between 
normal and stable behavior occurs at a = 2. Further, from (2.5) for z = A 2 , 
we conclude that E[Zi] = 1 for Z\ = V* 2 + ... + V^ 2 . It follows from (2.6) 
that m(s) < 1 for all s > a and it can be checked that E[|Zi| 2 ] < oo. Thus, 
the sufficient condition of Proposition 1.1(c) applies and Z n —>• Z a.s. and in 
C 2 for a complex random variable Z with P(Z ^ 0) > 0. 
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By Theorem 1.2, the set of solutions to (2.3) is given by all laws of random 
variables of the form 

Y w + aZ (2.7) 

where (Y t )t>o is a strictly (U, a)-stable Levy process independent of (W,Z). 
Here, U = {e X2t : t £ R}. For b = 27, this group is depicted in Figure 1.2a. 
The solution of interest can be singled out by moment properties using Remark 
1.4 and is X = ^iZ. This gives in particular a positive answer to the question 
posed in [29, Remark 4.5] about the existence of further solutions with infinite 
second moment. 

2.1.2 Cyclic Polya urns 

Consider an urn containing finitely many balls of b different types, 1,. .., b. 
At each step, a ball is drawn and placed back into the urn together with an 
assortment of new balls, the types of which depend on the type of the ball 
drawn. Such a scheme is called a generalized Polya urn. If the replacement 
rule is such that if a ball of type k is drawn, then it is placed back into the 
urn together with a ball of type k + 1 if k < b and of type 1 if k = b, the urn 
is called cyclic. 

Let R n} k be the number of balls of type 1 in a cyclic urn after n steps 
when starting with exactly one ball of type k and no other ball. We have 
E[i? n ,fc] = t + 0(1) as n —>■ oo, see e.g. [53, Lemma 6.7]. 

If b < 6, then R n> /- — j, suitably scaled, is asymptotically normal. If b > 7, 
let £ = exp(27ri/6) = £ + vq be a primitive 6th root of unity. It has been 
shown with martingale methods [46,69] and via the contraction method [53, 
Section 6.3] that n~^(R ni k — f) has an asymptotic periodic behavior (similar 
to Eq. (2.2)) that is governed by the law of a random variable X which is 
the unique non-degenerate solution with expectation 2 /(br(C, + 1)) and finite 
second moment of the equation 

X *= C/ c Xi + C (1 - ufx 2 (2.8) 

where Xi,X2 are i.i.d. copies of X that are independent of U which has the 
uniform distribution on [0,1]. 

Therefore, 


m(s)=E[\Ut\ s + \<;(l-U)t\ s ] = T ^- s , 

see Figure 1.1b for a plot of s e> m(s). Thus a = l/£ and a £ (1,2) iff 
£ = cos(27 x/b) > 1 iff b > 7. In particular, the phase transition at b = 6 is 
visible here as the phase transition between normal and stable behavior occurs 
at a = 2. 

(Al)-(A3) are readily checked to be valid in the present context. Further, 
E [Zi] = E[f7 c + C(1 - H) c ] = (1 + C)E[17 C ] = 1. 
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Since m(2) < 1 and E[|i?i| 2 ] < 4 < oo, the sufficient condition in Proposition 
1.1(c) is fulfilled and we conclude that Z n ^ Z a.s. and in C 2 as n —> oo for a 
random variable Z with E [Z] = 1 and E[|Z| 2 ] < oo. We have U = {£ fc e^ : 0 < 
k < 6, t £ R}; Figure 1.2b is a depiction of U in the case b = 7. It follows that 
the whole spectrum of solutions given in (1.13) appears. The special solution 
X appearing in the description of the limiting behavior of R Ut k — is the 
unique solution to (2.8) with mean 2 /(br(C, + 1)) and finite variance. Since 
E[|1V| 2 ] = oo for any non-trivial (U, a)-stable Levy process (Y t ) t >o by Remark 
1.4, it is X = 2Z/{br(C+ 1)). 

2.1.3 Asymptotic size of fragmentation trees 

In Kolmogorov’s conservative fragmentation model [15,54] an object of mass 
x = 1. say, is split into b parts with respective masses 0 < Vl, ..., 14 < 1 where 
b > 2 is a fixed integer and 14,...,14 are random variables with 14 +... + 14 = 
1 a.s. The splitting procedure is repeated with the resulting objects using 
independent copies of the splitting vector ( 14 ,..., 14 ) to determine the relative 
sizes of the emerging objects. Janson and Neininger [47] investigated the size 
N(e) of the random fragmentation tree the vertices of which correspond to all 
objects created in the fragmentation process that have mass strictly > e for 
some given e > 0. They showed that the asymptotics of N(e) are coded in the 
function if that maps z G C to E[J ^ =1 Vc] (whenever the expectation exists). 
To be more precise, denote by 1 = Ai,A 2 ,A 3 ,... the roots of the equation 
if{z) = 1 with the convention that 1 = Re(Ai) > Re(A 2 ) > Re(Aa) > .... 
Then, under suitable assumptions, when Re(A 2 ) < |, N(e) suitably shifted 
and scaled, converges in distribution to a centered normal. On the other hand, 
when Re(A 2 ) > N(e) exhibits a periodic limiting behavior governed by a 
complex-valued random variable X , with finite second moment and a fixed 
expectation 7 € C, satisfying the distributional equation 

X '= (2.9) 

3 =1 

where X 1 ,... ,Xb are i.i.d. copies of X and independent of (Vj,..., 14). 

This equation is in the scope of our analysis. Indeed, letting Tj := V) Az for 
j = 1 ,... ,b, and Tj = 0 for j > b , we have 

r -1 r b 1 

m(s) = E ^|T,| S = E J2v* e(X2)s = V’(Re(A 2 )s). (2.10) 

'-j> 1 2 '-1=1 

I 11 particular, m(a) = 1 iff a — l/Re(A 2 ). Again, the phase transition between 
normal and stable fluctuations is reflected in the equation since a < 2 iff 
Re(A 2 ) > \. Further, the conditions (A1)-(A3) are easily checked to hold. 
Since m( 2) < 1 and 

62-1 r b 

E y / 2 < bE ^y/ Re(A 2) < 1, 

L j=1 J l j=1 J 
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the sufficient condition of Proposition 1.1(c) is fulfilled and hence the mar¬ 
tingale (Z n ) nG f$ 0 converges a.s. and in C 2 and the special solution X used to 
describe the limiting behavior of N(e) is 7 Z. The general form of solutions is 
given by (1.13). 

It is assumed in [47] that each Vj has an absolutely continuous component, 
hence we are in the continuous case; in fact, U = {e X2t : t G R}. If the values 
of A 2 are the same, we obtain the same class of (U, a)-stable Levy processes as 
in the case of 6 -ary search trees, but the law of W will depend on the explicit 
distribution of Vi,..., 14, not only on their support. 


2.1.4 Biggins' martingale with complex parameter 
and complex Gaussian multiplicative chaos 


Eq. (1.10) arises naturally in the context of branching random walks: Consider 
an initial ancestor at the origin with children placed on R according to a point 
process Z on R with E[£(R)] > 1 (supercritical case). Each child produces 
offspring with positions relative to its location given by an independent copy 
of Z, and so on. Denoting by (<S(i;))|„|_ n the positions of the nth generation 
particles, consider the Laplace transform with complex parameter A of the 
random point measure formed by the nth generation particles, 

M n (X) ■■= J2 ■ 

\v\=n 


If m(A) := EEh=i e AlS ^] is finite, then E[M„(A)] = m(A)", and 


W„(A) 


Mnj A) 
m(A) n 


is a complex-valued martingale with E[Wi(A)] = 1, called Biggins ’ martingale, 
see [17]. Sufficient conditions for the convergence of W„(A) to a nondegenerate 
limit W(A) are studied in [18]. Upon defining Tj := e~ xs ^ /m(X) and using 
the same shift notation as for the weighted branching process, one obtains that 

. , e -AS(j) _ , 

W) = E^7TT [W(A) ]i = E T ^ W ( A )]i a ' s - 

j> 1 1 ’ 3> 1 


Thus, W(A) is a solution to (1.10), in particular, W„(A) = Z n in our notation. 

The sufficient conditions for the /^-convergence of W n (A) from [18, The¬ 
orem 1] translate as follows: (2.1) there is equivalent to E [(^) J . >1 |T)|) 7 ] < 00 
for some 7 G (1,2], while (2.2) equals m((3) < 1 for some /3 G ( 1 , 7 ]. These 
imply the sufficient conditions of Proposition 1.1(c). 

It is an important open problem to find equivalent conditions for the con¬ 
vergence of Z n to a nondegenerate limit, for it may also provide educated 
guesses in the theory of complex Gaussian multiplicative chaos. This is the 
complex analogue of real Gaussian multiplicative chaos, which was introduced 
by Kahane [49], see also [70] for a recent review and more details. Complex 
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Gaussian multiplicative chaos is a complex random measure A4 7 " 8 (with pa¬ 
rameters /3 ,7 > 0) on R, which is obtained via a limiting procedure from 
regularized measures AdJ’P. The question is about the correct renormalization 
needed to obtain convergence. Three phases (Phases I, II and III) appear, see 
Figure 1 in [56]. The suitable scaling can be guessed from the behavior of Big¬ 
gins’ martingale with complex parameter. A particular instance, which was 
studied by Madaule, Rhodes and Vargas [61] in order to provide intuition for 
the behavior on the boundary between phases I/II (7 £ (1/2,1), 7 + /3 = 1), is 

M n ( 7,(0) := ^ exp ( - 7>S(u) + i/3\/21n2 t S'(u)), 

\v\=n 


where (<S(u))„ € v, {S'{v)) v ^y are independent branching random walks with 
binary branching and i.i.d. displacements with normal laws with mean 2 log 2 
and variance 2 log 2 for S{v) resp. mean 0 and variance 1 for S'(v). 

As described above, if W n '■= (E[AI„]) _1 AI n converges to a limit W, then 
this limit is a solution to a smoothing equation, and equal to the particu¬ 
lar solution Z. In the setting of [61], one has E[A4 i(7 ,/3)] = 1 and m(s) = 
exp ((sj— l) 2 log 2), hence a = I/7 £ (1,2). The phase transition is reflected in 
the fact that m'(a) = 0 on the boundary between phases I/II. In this case, the 
sufficient conditions of Proposition 1.1 do not hold, nevertheless it is proved 
in [61, Theorem 1] that W n converges to a nontrivial limit. 


2.2 Application of Theorem 1.5 

We end this section with the study of the example described in Section 1.2.2. 
In contrast to the examples above, the relevant solution will be given by Yw ■ 
Bassetti and Matthes [14, Section 6.2] study the equation 

V '= LV\ + RV 2 , (2.11) 

where V, V), V/ are i.i.d. random vectors in R 3 , independent of the random 
pair (L, R) of similarities which satisfies 

m( 2) = E[ ||T|| 2 + ||I?|| 2 ] = 1 and m(p) = E[ \\L\\ P + ||i?f ] < 1 

for some p £ (2,3). Thus (A1)-(A3) are satisfied with a = 2, and all solutions 
to (2.11) are given by Theorem 1.5. Since the equation is homogeneous, W* 
vanishes, while Y\ can be any centered multivariate normal random variable 
with a covariance matrix that is invariant under conjugation by elements of 
O. 

Bassetti and Matthes further assume that L = IA for independent random 
variables l £ R and A £ 0(d), and that the law of A T x dominates the volume 
measure on § d_1 for every x £ § d_1 . This implies that O acts transitively 
on S d_1 and hence that scalar multiples of S = are the only possible 
choices for the covariance matrix of Y\. As Proposition 4.4 below shows, the 
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weaker assumption that there is no U-invariant proper subspace readily implies 

S = I d . 

Finally, Z = Z w ^ 0 if and only if Lw + Rw = w a.s., which corresponds 
to the conservation of momentum. In [14], only centered solutions were con¬ 
sidered, the physical interpretation of which is that the centre of gravity does 
not move. Here, we see that the centre of gravity may have drift Z w , which 
corresponds to the validity of Newtons first law in this context. 

In recent papers by Dolera and Regazzini [34] and Bassetti et. al. [13], 
solutions of the Boltzmann equation for Maxwellian molecules in R 3 have been 
studied directly, rather than its simplifications like the Kac caricature. These 
steady state solutions cannot be written directly as solutions to smoothing 
equations, but the techniques employed seem to be very similar. We hope 
that our results allow for a better understanding; in particular under which 
conditions rotation invariant solutions appear. 


3 Proofs of the main results 

In this section, we prove our main results: Proposition 1.1, Theorem 1.5 and 
Proposition 1.6. At the beginning, we collect the relevant notation and intro¬ 
duce tools and concepts which are used in the proofs below. 


3.1 Notation 

3.1.1 Vector spaces, sets, matrices, etc 

We work in the e?-dinrensional space R d . We think of an element x £ R rf as 
a column vector. We write x T for the corresponding row vector, ei,..., 
denote the canonical basis vectors of R d . By (•,•), we denote the standard 
Euclidean scalar product on R d , that is, ( x,y) = x T y for x, y £ R d . We write 
| a: | for \J (x, x), the Euclidean norm of x. For a set B C R d , dB denotes 
the boundary of B, B - 1 = {y £ R d : (x,y) = 0 for all x £ B} denotes the 
orthogonal complement of B in R d . B r := {a: £ R d : |ar| < r} denotes the 
ball of radius r centered around the origin, r > 0. For a given real d x d 
matrix A , we write Ajj for the coefficient in the ith row and the jtii column 
of A and ||A|| := sup| x i =1 |Aa;| for its norm. A T ,tr(A) and det(A) denote the 
transpose, the trace and the determinant of A, respectively. For A £ R, we 
set E\(A) = {x £ R d : Ax = Aa:}. We write Ik for the k x k identity matrix, 
k £ N. 


3.1.2 Probability spaces, expectations, etc 

Throughout the paper, we fix a probability space (12, A, P) which is large 
enough to carry all random variables appearing in the paper. By E[-] and Var[-] 
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we denote expectation resp. variance with respect to IP. We also consider expec¬ 
tations of random vectors and random matrices, which are defined componen¬ 
twise. For a random vector Y, we write Cov[Y] := E[(Y — E[Y])(Y T — E[Y] T )] 
for the covariance matrix of Y with respect to P. denotes the indicator 
function of a set A and we write E[Y; A] for E[Y 1^]. Further, Cov[Y; A] is the 
covariance matrix of the random vector Y 1a- 


3.1.3 Relevant groups 

The following groups are of relevance. S (d) C GL(d, R) is the group of similar¬ 
ity matrix, i.e., scalar multiples of orthogonal matrices. U denotes the smallest 
closed subgroup of S (d) that covers {Tj : j = 1 with probability 

one. In the complex case, we identify U with a subgroup of the multiplicative 
group C* of C. Similarly, O is the smallest closed subgroup of O (d) which 
contains the random set {Tj/ ||T)|| : j = 1 with probability one. In 

the complex case, we identify O with the smallest multiplicative subgroup of 
§ = {zGC:|z|=:l} which contains the random set {Tj/\Tj\ : j = 1,..., N} 
with probability one. In this case, either O = {e 27nfc / m : k = 0 ,... ,m — 1} 
for some m £ N or O = S. Finally, let G be the smallest closed multiplica¬ 
tive subgroup of R> that covers the random set {||T)|| : j = 1,... ,N} with 
probability one. Equivalently, G is the image of the group U under the homo¬ 
morphism u i-)- Hull- There are three possibilities: (C) G = R>; (G) G = r z for 
some r > 1; (T) G = {1}. The trivial case (T) is excluded by (A3). We refer 
to (G) as the geometric or r-geometric case and to (C) as the continuous or 
non-geometric case. 

More information about the structure of U is provided in Section 4.1. In 
particular, U = A \j k Cu for a one-paranreter group An which is isomorphic 
to G, and Cu =UnO(d). Note that in general, Cu C O. 

Throughout the paper, we call a measure g on R d G-(left)-invariant for 
a group G of matrices if g{g~ 1 B) = g{B) for all g G G and all Borel sets 
B C R d . 


3.1./ Weighted branching 

We set T n := cr((C(v),T(v)) : |i>| < n) and T := a(T n : n £ No). (C,T) 
serves as an abbreviation for the family ((C(v), T(v))) ve y. We further assume 
that on the basic probability space, a family X := (X v ) ve y of i.i.d. random 
variables is defined which is independent of T. We do not specify the law of 
X := X 0 here, but typically, X will be a solution to (1.1) or (1.7). 


3.2 Characteristic functions 

It is natural to approach Equation (1.1) via characteristic functions. Indeed, 
the distributional equation (1.1) for an Revalued random variable X is equiv- 
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alent to the functional equation 


<j>(x) = E 


e^ c> n 4>{T]x) 

3> 1 


x€R d 


(3.1) 


for the characteristic function <f>{x) = Efe'^A)] of X. In the homogeneous 
case, the functional equation takes the simpler form 


4>{x) 


E 


W^Tjx) 


L j> i 


x G R d . 


(3.2) 


Solving (1.1) is equivalent to finding all characteristic functions <j> of R d -valued 
random variables which satisfy (3.1). 


3.3 Proof of Theorem 1.5: The direct inclusion 

We are now ready to prove the direct inclusion of Theorem 1.5. Here, with 
the direct inclusion, we mean the assertion that any distribution of a random 
variable of the form (1.16) is a solution to (1.1). 


Proof (of the direct inclusion of Theorem 1.5) Consider the situation of The¬ 
orem 1.5 and let X = W* + Yw + Z where (Y t ) t > o is a strictly (U, a)-stable 
Levy process independent of the family (C, T), in particular independent of 
(W*,W,Z). Denote the characteristic function of X by <j> and the character¬ 
istic exponent of Y\ by T, that is, Efe 1 ^’ 1 *)] = exp (t\P(x)). By assumption, T 
satisfies (1.12) with b(u) = 0 for u G U. Using the independence of (IT*, IT, Z) 
and {Yt)t> o, we conclude that 


4>(x) = E exp (i(a;, W*) + i(cc, Z) + W\P(x ))], x G (3.3) 


As ([W*]j, \W]j, [ Z]j ) is a copy of (W*,W, Z), (3.3) still holds when (W*,W, Z) 
is replaced by ([W*]j, [W]j, [Z]j), j G N. Furthermore, since ([W*]j, \W]j, [. Z]j ) 
is independent of T\. we conclude that 


cf>(Tjx) = E[exp(i(T/*, [W*} 3 ) + i (Tj x, [Z} 3 ) + [W]^{Tj. x)) \ X,} 

= E [ exp(i(x, Tj \W*\ 3 ) + i (x,T 0 [Z\ 0 ) + \\Tjf [W} 0 T{x)) \ X,} 
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a.s. for every x £ R d , where in the last step we have used (1.12) and the fact 
that Tj £ U a.s. Consequently, for every x £ R d , we infer 


e^n <KTjx) 

3> 1 J 



= E 

e i{x ’ G) E 

L j>i 

r e i<x,T i [W*] 3 >+i<x,T 3 [Z]j) + ||T J |rhy] ; ,-f(x) 

X 1 

= E 

[ e i(x,c )+ E,>i 

(i(x,T J [lV*]j)+i<x,T J [Z]j> + ||T 3 |nw]^(x))' 


= E 

Li^.C+E^i hriil+ib.EjM T J [Z] j )+j: jS:1 \\T j \\ a [W] J *(x)'] 

= E 

r e i<x,W*)+i<x,Z)+rV0'(a:)l _ 



i.e., 4> solves (3.1). Hence, X is a solution to (1.1). □ 

Most of the remainder of this paper is devoted to the proof of the converse 
inclusion of Theorem 1.5 and the description of the class of strictly (U, ca¬ 
stable Levy processes. We begin with a short section on a common technique 
in the theory of branching processes, an exponential change of measure. 


3.4 Exponential change of measure 


Recall from Section 1.3.3 the definition of the weighted branching process. 
We define the associated random walk (L n ) ne n 0 on S (d) by the many-to-one 
formula 


E[f(L 0 ,...,L n )\ 


E 


E 

\v\=n 


\\L(v)\\ a f((L(v\ k )) k =o,...,n) 


(3.4) 


for all nonnegative Borel-measurable functions / : S (d) n+1 —> R>. (A2) im¬ 
plies that the law of (L n ) ne n 0 is a proper probability measure. From this 
definition, it can be checked that (L n )neN 0 is a multiplicative random walk 
on the group U C S(d). We define S n := — log||L„|| and O n := L n /\\L n \\. 
(Sn) n ^n 0 is a standard random walk on R, while (O n ) n eN 0 is a multiplicative 
random walk on O (d). The step distribution of (S' rl ) n€ N 0 is given by 


nsi £ •) 


r N 


E 


Ell T jl| Q<5 - lo gll T hl(‘) 


3=1 


(3.5) 


Consequently, when (A3) holds, E[Si] = — m'{a) £ (0, oo). 

Later on, we will use that (3.4) remains valid under certain stopping rules: 
Considering r(f) := inf{n £ No : S n > t and S k < t for all k < n}, (3.4) gives 


E[/(Ao, • • • ,L n )t { r (t)= n }] 
= E 


'y ' l|A(u)|| /((L('c|fc))fc=0,...,n) l{S(x)>i>S(x| fe ) Vfc<n} 

\v\=n 
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Summing over all n £ No and defining the coming generation at time t 
C(t) = {»£?: ||L(u)|| > 0 and S(v) > t > 5(v|fc) for all k < |u|}, 
we infer that in particular 


E[/(S r(t) _ L ,S rW ,O rW )] =E 




> 0 , 

(3.6) 


(3.7) 


L t >ec(t) 


for all nonnegative Borel-measurable functions / : K 2 x O (d) —> K>. See [55] 
for more information on stopping lines and further references. 


3.5 Multiplicative martingales 

Let X be a solution to (1.1) and denote the characteristic function of X by (/>. 
We will show that (f> is the characteristic function of a random variable of the 
form (1.16). As in previous works on fixed points of smoothing transformations 
[2,4,5,20,21,43], we make use of multiplicative martingales. As this technique 
is well-known by now, we keep the presentation short here and refer to the 
above references for more detailed expositions. For x £ R d , define 

M n (x) := exp(i(a;, W*)) ■ (f(L(v) T x), n £ N 0 . (3.8) 

\v\=n 

The fact that <j) solves (3.1) implies that (M„(a;))„ e N 0 is a complex-valued 
martingale. Since it is bounded by 1 in absolute value, it converges a.s. and in 
mean. We denote its a.s. limit by M{x) and note that 

(t>{x) = E [M(x)\, x£R d . (3.9) 

In order to determine <p(x), it thus suffices to determine M{x). 

We begin with a key lemma. We will sometimes write M(x) = M(oj,x) 
in order to make more transparent, whether we consider M(x) as a random 
variable, or x i —> M(ui, x) as a function of x. 

Lemma 3.1 There is an T-measurable set N with P(N) = 0 such that, for 
(jo £ N c , M n (u, ■) — > pointwise on Further, x i —> M(ui,x) is the 

characteristic function of a probability distribution on M. d for all oj £ N c . 

An obvious modification of the proof of Theorem 1 in [26] yields the result. 
We refrain from giving any details. 

Solving the inhomogeneous equation (1.1) can be reduced to solving the 
associated homogeneous equation (1.7) along the lines of [5, Section 5], a sketch 
of the reduction argument will be given in Section 3.11 below. For now, we 
restrict our attention to the homogeneous case and assume that (7 = 0 a.s. 
Then (3.8) takes the simpler form 

M n (x) = 4 >(L{v) t x), n € N 0 . 

\v\=n 


(3.10) 




26 


M. Meiners and S. Mentemeier 


We claim that for all ui from a set of probability one, M (oj. •) is the characteris¬ 
tic function of an infinitely divisible law. Since sup|„i_ n ||L(u)|| —> 0 a.s., we can 
assume without loss of generality that the set N from Lemma 3.1 is such that 
N c C {sup|„| =n ||L(u)|| —> 0}. Now pick an w £ N c . We view M n (u,-) as the 
characteristic function of 'Yl\ v \— n L{v)X v conditional given (L(v)) ve y. Thus, 
M(u), •) is the limit of characteristic functions of the row sums in a triangular 
array which is independent and infinitesimal. Hence, it is the characteristic 
function of an infinitely divisible law and we can write M(x) = exp (S'(a;)) for 
a random characteristic exponent \P(x ) which is of the form 

*{x) = i (W,x)-'^+ J (e^-l-i {x,y)l [0!l] (\y\)) u(dy) (3.11) 

for W' = W'{ T) £ R d , a covariance matrix S = S(T) and a Levy measure 
v = v(T) on R d (see [52, p. 290]). 

That W',2J and v are indeed functions of T comes from the following 
representation, valid for triangular arrays (see [52, Chapter 15]): Using the 
convention / {h<|x| < 1} = “ / { i<|*|<h } when h > ^ d holds on 

W' = W h + [ xv{dx) (3.12) 

J {/l<|a;|<l} 

where W h is defined by 

W h := lim Y E[L(v)X v :\L(v)X v \ < h \ T] (3.13) 

n—>oo • ^ 

\v\—n 

for every h > 0 with v({\x\ = h}) = 0, 

S = 'E h - [ xx T v{dx) (3.14) 

J {h<|aj|<l} 


where is dehned by 

S' 1 := lim y Cov\L(v)X v ; \L(v)X v \ <h\T], (3.15) 

n—>oo ' ^ 

\v\=n 

and 

[ f(x)v(dx) = lim y [ f{L(v)x) F(dx) (3.16) 

J n—t oo • ^ J 

|t>|=n 

for all continuous functions / with compact support on R d \ {0} (R d denotes 
the one-point compactification of R d ). Eq. (3.16) also yields that v is a ran¬ 
dom measure (see [51, Lemma 4.1]), i.e., the mapping T f f(x)v(T, dx) is 
measurable for every nonnegative Borel measurable function / on R d \ {0}. 
In order to show that W’ and S are random variables as well, we need some 
more preparation (the problem is to choose ft in a measurable way). 
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Right now, we can use that W'. E and v are functions of T and apply the 
shift operator [•]„. Arguing as in the proof of Lemma 4.3 in [5], we conclude 
that on N c 

M(x) = [M] v (L(v) t x) for all x G R d and n G No- (3-17) 

| V | —71 

Now using (3.17) in (3.11), we conclude that for all n G No, on N c 

m = iE wum T x) - E EtEmEEi 

|t>|=n |u|=n 

+ E /( ei<i(l,)Tx,!/> - 1 - i<i(«) T *,|/)l [ 0 il] (|j/|)) M„(d y) 

\v\=n 

= i( E L i v )[ w ']v, x ) ~\ x * E L{v)[E\ v L(v) T x (3.18) 

|i;|=n \v\=n 

+ E / (e l{xMv)v) - 1 - i{x,L(v)y)l [0tl] (\L(v)y\)j[i/] v {dy) 

\v\=n 

- E y , ( i <*> L ( u )2/) ]l [o,i](bl)- i < a; »- z: '( w )y) 1 [o > i](l-C'(%l))[i / ]«( d J/)- 

\v\=n 

Since the last term contributes to the random shift, using the uniqueness of 
the Levy triplet, we get: 

J f(y)v(dy) = j f( L ( v )y)W\v(dy) on N c (3.19) 

l^l^n 

and E = ^ L(v)(E] v L(v) J on N c (3.20) 

\v\=n 


for all n G No and all nonnegative Borel-measurable functions / on \ {0}. 
We will use (3.19) and (3.20) to determine v and E, respectively. 

The next lemma gives an important estimate for v and will allow us to 
infer the measurability of W' and £. 

Lemma 3.2 There is a multiplicatively r-periodic (if & = r z ) or constant (if 
G = K >y ) function i) : [0, oo) — > (0, oo) such that t H > t)(t)t a is nondecreasing 
and such that on N c , 

v(B\ x j) = Wfi(|a:|^ 1 )|a:|-“ for all x G R d \ {0}. (3.21) 

Proof Define / : [0, oo) —> [0,1] by 


E[exp 

1 




for t > 0, 
for t = 0. 


/(*) ■ 
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f is decreasing in t and continuous at 0 since i 4- 0 as t 4- 0. The most 
important property of /, however, is that it solves the functional equation of 
the smoothing transformation as studied in [2]. Indeed, using (3.19) and the 
fact that L[y) = ||L(u)|| 0(v) for an orthogonal matrix 0(v) and 0(v)~ 1 B = 
B$ for all r > 0, we get 


/(*) = E 

= E 


ex P (- l u ]v( L ( v ) ls t-0 

' \v\=n 


n ex p( 

|p|=n 




= E 


n fmv)\\t) 

| v\=n 


Consider the limit Mj of the multiplicative martingale associated with /, i.e., 
Mf(t) = limn^oo ]~[|u|=n f(\\L{ v )\\ t), t > 0. By Theorem 8.3 in [2], M f (t) = 
exp(— Wtf(t)t a ) a.s. for all t > 0 and some function f) with properties as above. 
By the arguments given in the proof of [5, Lemma 4.8], 

v{B \ x |) = W^x]- 1 )^]- 01 for all x £ R d \ {0} on N C .D (3.22) 

Since (3.21) holds, we can pick some 0 < h < 1 such that fi is continuous 
at h, in particular, v(dBh) = 0 on N c . Using this h in (3.12), (3.13), (3.14) 
and (3.15) implies the asserted measurability statements for W' and S. 

We call a Levy triplet (W', S, v) with these measurability properties an 
^-measurable Levy triplet. We summarize the results of the above discussion 
in the following proposition. 


Proposition 3.3 Assume that (Al) and (A2) hold. Let X be a solution to 
(1.7) and denote its characteristic function and distribution function by <f> and 
F, respectively. Denote by (M. n (x)) n ^ 0 , x £ the multiplicative martingales 
associated with X and defined by (3.10). Then, on an F-measurable set N c 
with P(AT) = 0, 

M n (x)—>M(x) as n —> oo (3.23) 

for all x £ R d . For u> £ N c , x i—>• M(cj,x) is a characteristic function and 
possesses a representation M(x) = exp(F(x)) for all x £ R d where T> is given 
by (3.11). {W', S, v) is an F-measurable random Levy triplet satisfying (3.12)- 
(3.16). 


Next, we state some consequences of Proposition 3.3 concerning the tail 
behavior of solutions to (1.7). These will be useful when determining v and S. 


3.6 Tail estimates 

If X is a solution to (1.7), then 

limsupt Q P(|X| > t) < oo. 

t->-oo 


(3.24) 
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If, additionally, the random Levy measure v of the limit of the multiplicative 
martingale M associated with X (or its characteristic function <f) vanishes 
a.s., then the stronger estimate 

limsupf Q P(|X| > t) = 0 (3.25) 

t->-oo 

holds. The derivation of these estimates can be carried out along the lines of 
[43, Lemma 4.7] and [5, Lemma 4.9], we refrain from giving more details here. 
As a consequence of (3.24), we obtain the following inequality for 

Lp(t):=E[\Xf-,\X\<t] (3.26) 


with (3 > a: 


limsupf Q ^Lp(t) < limsupf a & f (3x^ 1 P(|X| > x) da: (3.27) 

t—>oo £—>oo Jo 

< C lim sup f“~' 9 [ x^-'dxKoo 

t —>oo Jo 

where C is a constant depending on a, /3 and the law of X (via (3.24)). In the 
case where the random Levy measure v vanishes a.s., using (3.25) instead of 
(3.24) produces the stronger estimate 

lim t a -^L p {t) = 0. (3.28) 


3.7 Determining v 

The characterization of v is given by the following lemma. 

Lemma 3.4 Assume that (A1)-(A3) are in force and let v be a F-measurable 
random Levy measure. Then v satisfies (3.19) if and only if v = Wv a.s. for 
a deterministic (U, a) -invariant Levy measure D, i.e., satisfying 

D(gB) = \\g\r a u(B) (3.29) 

for all g € U and all Borel sets B C R d \ {0}. Further, v = 0 if a > 2. 

Proof We first prove the sufficiency. To this end, suppose v = Wv a.s. for a 
deterministic Levy measure v satisfying (3.29). It suffices to check validity of 
(3.19) when / is the indicator function of an arbitrary Borel set B C R d \ {0}. 
Using (3.29) and then (1.9), we obtain 

E WUm^B) = E [W} v h(L(v)- 1 B) 

|l)| =n |u|=n 

= E 11 ^) 11 “ [W] v u(B) = WD(B) = v(B) 

\v\=n 


a.s. 
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For the converse implication, assume that v is an ^-measurable random 
Levy measure satisfying (3.19). 

The idea of the proof is as follows: We define v := E[i/] and prove that 
9 ^ llffl| -Q v(g~ 1 B) is a constant function in g (this necessitates proving that 
{ui,g) i —y v(g~ 1 B) is product-measurable). Then we show that the identity 
v{B) = WD(B) a.s. holds simultaneously for enough sets B (namely, a count¬ 
able generator of the Borel cr-field which is closed under finite intersections) 
on a common P-null set. 


Step 1: We start by proving the measurability statement: For every Borel 
set B C \ {0}, the mapping (u,g) i-)- v{g~ 1 B) from fi x § (d) to R> 
is T ® $ 8(S(d))-measurable where ®(§(d)) denotes the Borel cr-field on §(d). 
Let B C R d \{0} be closed, B l / k := {y £ R d : \y—z\ < f for some z £ B} and 
fk be a continuous function satisfying 0 < fk < 1, /a-|b = 1 and /fc|( B i/ky = 0 
(such a function exists, for instance, by Urysohn’s lennna). Then fk —> 1b as 
k —> oo. By (3.16), 

[ fk(gx)u(dx) = lira V [ f k (L(v)gx)F(dx) 

J n—>•oo z ' J 

|v|—n 

By Fubini’s theorem, the right-hand side is T (g> *B(§(d))-measurable and by 
the dominated convergence theorem, the left-hand side tends to v{g~ k B) as 
k oo. Hence (w,g) >->• v{g~ 1 B) is T <g> *B(§(d))-measurable. For every fixed 
r > 0, this extends to the Dynkin system generated by all closed sets B C B x , 
hence to all Borel measurable sets B C \ {0}. 


Step 2: Now we introduce the sets that will form the countable generator. We 
define, for x £ R d and e > 0, 


II = { y £ R d : \y\ > |x|, 


V_ _ x_ 

\y\ M 


< e 


We will determine I'C^/ixp)- First notice that, for o £ ©(d), 


oy x 


-if v x \ 


JL _ -i JfL 

\oy\ \x\ 


V \y\ M ) 


\y\ ° M 


and, therefore, for a similarity g = ||g|| o with ||g|| > 0, 

s~'n = e R,i: ,!, i a ]jjj| • | jfi - 0-1 pj 

where in the last step, we have used that g _1 = 
this in (3.19) with / = 1/e gives 

x/\x \ 2 


- 6 c Ig~ l x I g T x/\\g \\ 2 


I —1 „T 


(3.30) 

° T = 9 T / \\g\\ 2 - Using 


u (Ix/\x\ 2 ) — X] M«(U(^) 1 Il/\x\ 2 ) — X] Hr'(-^L(v) T x/|L(u) T x| 2 ) a - S - (3-31) 

|t»|=n |u|=n 
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Step 3: For x E \ {0}, define <P e (x) := v(Il^ x ^)\x\ a . By (3.21), 1%^ C 
B \ x |-i implies 

K 7 x/|*p) < H B \ x \-i) = Wt)(\x\)\x\ a a.s. (3.32) 

Consequently, with h* := sup t>0 f j(t), we have 

sup & £ (x) < WA* a.s. (3.33) 

xeR d \{0} 


Further, by (3.31), 

<F e (a:) = E ||L(u)|| q [ | F e ]^(L(u) T a;) for all x £ R d \ {0} a.s. (3.34) 

\v\=n 

For x E M d \ {0}, define ip €tX : U — > R> via 

A,x(g) = R{&e{g T x)}, g & U. 

is measurable (since (ui,g) H > v(g~ l B) is product-measurable), nonnega¬ 
tive and bounded (due to (3.33) and the fact that E[W] = 1). Further, as a 
consequence of (3.34) and the change of measure (3.4), we have 


A,x(d) = K [^e(g T x)} = E 


= E 


E \\L(v)\\ a [* e ] v (L{v) T g T x) 

|u|=n 


E \\ L ( v )\\ a ^x{ 9 L ( v )) 

|u|=n 


= E[lp^ x {gL n )] 


for all g E U. In other words, ip e , x solves a Choquet-Deny functional equation 
on the (possibly non-Abelian) topological group U. From the Choquet-Deny 
lemma (see Lemma A.l in the appendix) we infer that ip e<x is P (L. n E -)-a.s. 
constant on U. Consequently, there is a constant c £>x > 0 such that, with 
U = {ip e ,x = c e,x}, P (L n E U) = 1 for all n E Nq. Notice that this implies 


E 


E II L(v)rML(v)) 


>l=" 


P (L n EU) = 1 


and hence P (L(v) E U for all |t>| = n with ||L(w)|| > 0) = 1. This together 
with the martingale convergence theorem and (3.34) yields 

#«(*)= lim E[<F e (a;)|J r n ] = lim E 11^)11“ E[[^ e ] v (L(v) T x)\JF n ] 

n —>oo n —>oo L ' 

\v\=n 

= lim E \\L(v)\r^, x (L(v)) = Wc e , x a.s. 

n —>oo z ' 

\v\=n 


In terms of , , 2 ), this reads 


K^/lzp) = Wc e , x \x\ a 


a.s. 
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Setting u(-) := E[i/(-)] and taking expectations in the above equation gives 
f>{I* x/ \ x? ) = c e>x |x| a and hence 

v Ax/\x\ 2 ) = Wv(Il/\x\ 2 ) a.s. (3.35) 

Step 4: Using (3.30), we conclude that, for every jeD, 

U (9~ lT x/ 1*1=0 = U ( I g T x/\g T x\ 2 ) = WAA9 T )\9 T x\ a = llfiff AA/\x\ 2 ) a - s - 

and hence, upon taking expectation, we arrive at the transformation formula 
(3.29) for P, valid for all g G U and sets JU. , 2 . (3.29) extends to finite in¬ 
tersections of the form I = 7^| , 2 D ... D /\ x p for x = (xi, ... ,x n ) G 
(K d \ {0}) n and e = (ei,...,e„) G R>. To be more precise, first notice that 
in such an intersection, one can assume that |xi| = ... = \x n \ (otherwise, let 
r := min” =1 \'Xj\ and notice that I does not change when Xj is replaced by 
rxj/\xj\, j = 1,... ,n). (3.32) then carries over with 7W , 2 replaced by 7 and 
|x| replaced by \x\\. The obvious counterparts for 7 of (3.30) and (3.31) hold 
true and one can define !F e (x) := v{I) |xi| -a . The remaining arguments apply 
almost without changes and give 

v(I) = Wv(I) a.s. (3.36) 

and the counterpart of (3.29). Now notice that the finite intersections of sets 
I x /\ x \2 with x G Q d \ {0} and rational e > 0 form a countable generator of the 
Borel (7-field on and on a set of probability one, v(I) = Wv{I) for all sets 
from this generator. It follows by an application of the uniqueness theorem 
from measure theory that v(B) = Wv{B) for all Borel sets B C \ {0}. 
Analogously, (3.29) extends to all Borel measurable B CM d \ {0}. 


Step 5: It remains to show that v is a Levy measure and that v = 0 if a > 2. 
For the former, notice that P is a measure on \ {0} and that 


J (|x| 2 A 1) v(dx) < oo iff 


J (|x| 2 A 1) v(dx) < oo a.s. 


by (3.36). For the proof of the second claim, notice that taking expectations in 
(3.21) gives D(B c t ) = 6 (t~ l )t~ a for allt > 0 for a multiplicatively G-periodic 
function ft such that b [t)t a is nondecreasing in t. Because of these properties 
we have ft* := inf 4> o b(t) = 0 iff ft = 0. Now assume that a > 2. We show that 
then ft* = 0. Indeed, since D is a Levy measure, integration by parts gives 


oo > 


v(&x) = f 2tv(Bt \B\)At 
J o 

= f 2tt)(t~ 1 )t~ a dt — ()(1) > 2ft* f t 1 ^ a dt — b(l) 
Jo Jo 


'{ M 2 <i> 


which implies that ft* =0. 


□ 
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Any Levy measure V satisfying the invariance property (3.29) can be fac¬ 
torized into a radial and spherical part, similar to the decomposition valid 
for Levy measures of stable laws. The details are given in Proposition 4.3 in 
Section 4 below. 


3.8 Determining £ 

Now we use the identity 

£ = L{v)['Z} v L(v) t a.s., (3.20) 

\v\=n 

to determine E. 

Lemma 3.5 Assume that (A1)-(A3) hold and let £ be a IF-measurable ran¬ 
dom covariance matrix. Then £ satisfies (3.20) if and only if £ = WE a.s. 
for a deterministic covariance matrix E satisfying 

E = oEo J (3.37) 

for all o € O. Further, E = 0 if a ^ 2. 

Proof Our first observation is that tr(£), the trace of £, is a nonnegative en¬ 
dogenous fixed point of the one-dimensional smoothing transform with weights 
(||T i || 2 ) i >i. To see this, notice that it follows from (3.20) and the fact that the 
trace is invariant under conjugations with orthogonal transformations that 

tr(£)=tr( £ L(v)[i:] v L(v)A = £ ||L(t;)|| 2 tr(0(w)[£]„0(t;) T ) 

= H- l <t)II 2 M s )]” a - s - 

\v\=n 

From Theorem 6.2 in [2], we conclude that tr(£) = cW for some constant 
c > 0 if a = 2, and tr(£) = 0 a.s., otherwise. 

Now recall that all eigenvalues of a covariance matrix, i.e., a positive semi- 
definite, symmetric matrix are nonnegative and that there is a basis of 
consisting only of eigenvectors of that matrix. Further, the trace of the matrix 
equals the sum of its eigenvalues. 

In the case a/2, this implies that all eigenvalues of £ equal 0 a.s., which 
means that the matrix itself vanishes a.s. This, in turn, proves the assertion 
in the case where a/2. 

In the case a = 2, we can conclude that tr(£) has finite expectation since 
E [IT] = 1. This implies that £^ is integrable for i,j = l,...,d since 

\Eij\ = |e,£ej| < ||£|| < tr(£) 

where for the last inequality, we have used that the norm of £ is the largest 
eigenvalue of £ while the trace of £ is the sum of its eigenvalues. Consequently, 
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£ := E[E] is a symmetric dxd matrix. We claim that £ = o£o T for all o £ O. 
To prove this, notice that since £ is symmetric and positive semi-definite, it 
has n < d distinct real eigenvalues Ai > ... > X n > 0 such that d\ + .. .+d n = d 
where dj denotes the dimension of the eigenspace E\. (A7) corresponding to the 
eigenvalue A j, j = 1,... , n. For arbitrary x £ E Ai(E), first apply (3.20) for 
n = 1 and then condition w.r.t. T\ to obtain 


Ai|a:| 2 


x T £x = E 


N 

^2\\T j \\ 2 x T 0 j £0]x 


i=i 


(3.38) 


On the other hand, x T Oj£Ojx < Ai |rr | 2 for all j = 1 a.s. since Ai 

is the largest eigenvalue of £. Using m( 2) = 1 and (3.38), we conclude that 
x T Oj£Ojx = Ai|x| 2 for all j = 1 a.s. Since x £ E\ 1 (£) was arbi¬ 

trary, we infer that, a.s., the restriction of Oj to E\ 1 (E) is an automorphism 
of E\ 1 (£), j = Using this fact, the above argument can be re¬ 

peated consecutively for the remaining eigenvalues of £ (in decreasing order) 
to conclude that, a.s., each Oj , j = 1,..., N is an automorphism of E \ k (£), 
k = 2,... ,n. Since every vector x £ R d can be written as a linear combina¬ 
tion of normed eigenvectors of the eigenvalues Ai,...,A n , we conclude that 
£x = Oj £ Oj x for j = 1,..., N a.s. and hence £ = Oj£Oj for j = 1,N 
a.s. Standard arguments then yield that £ = o£o T for all o £ O. Using this, 
the martingale convergence theorem, (3.20) and the convergence 1U„ —> W 
a.s., we conclude that 


E = lim E[E|j;] = lim V L(v)£L(v) J 

n —±00 n —±00 ' ^ 

\v\=n 

= lim y \\L(v)\\ 2 0{v)£0{v) t = W£ a.s. 

n—> 00 • ^ 

\v\=n 


□ 


3.9 The proof of Proposition 1.1 

In contrast to Eqs. (3.19) and (3.20), which allowed to determine v and £, 
there is in general no such equation for lU' due to additional contributions 
coming from the Levy measure, see Eq. (3.18). Therefore, we have to postpone 
the identification of W' until the next section, and will prove Proposition 1.1 
first. This is closely related to determining W': Indeed, suppose that Z is a 
solution to (1.10), in particular, Z is ^-measurable. Iteration of (1.10) yields 

Z = y L(v)[Z] v a.s. (3.39) 

\v\=n 


for all n 6 No- Then consider the multiplicative martingales (M„(cc)) rae N 0 
associated with Z (and its characteristic function <fi). As in the proof of [5, 
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Proposition 4.17], we infer from the martingale convergence theorem, 


M(x) = lint TT cf>(L(v) T x) 

n —>00 

\v\=n 


lim E 


exp 




= lim E 

n—too 


exp [ 1 ^ x, £ mvi, 

\v\=n 

= exp(i(a:, Z)) a.s. 


Tn 


lint E[exp(i(x, Z))\F n ] 

n—> 00 


Hence, Z = W' a.s. in (3.11) with <P denoting the characteristic exponent of 
M (in this case, £ and v vanish a.s.). Therefore, by Proposition 3.3, Z can be 
written in the form 


lim y E[L(v)[Z] v ;\L{v)[Z] v \<l\T n ] a.s. (3.40) 

n—^oo z ' 

\v\=n 


Further, (3.25) gives 


P(|Z| > t) = o(t a ) as t —> 00 . (3.41) 

This and the estimate (3.28) for Lp (see Eq. (3.26)) from Section 3.6 is every¬ 
thing we need to determine Z when a ^ 1. 

Proof (Proof of Proposition 1.1 in the case a ^ 1) (a) Let 0 < a < 1. Then, 
using (3.40), (3.28), (1.5) and W n W a.s., we infer 


\Z\ 


lim 

n—too 


E K[L(v)[Z] v -\L(v)[Z] v \ < 1 | F n \ 

\v\—n 


< limsup Y ||£(t;)||E[|[Z] l) |;|[Z] t ,| < IIL^H" 1 | P n \ 

n—¥oo . , 


= lim sup E ||L(w)|| Li(||L(u)|| 1 ) = 0 a.s. 

n —>00 . . 


(c) Suppose that 1 < a < 2, that Z satisfies (3.39) and that P(Z ^ 0) > 0. 
Notice that (3.41) implies that w := E [Z] is finite, moreover, Z £ C s for all 
s < a. By standard martingale theory and (3.39), for any 1 < s < a, 


Z = lim E[.Z|.7 r n ] = 


lim E 

n—> 00 


E 


\v\=n 


Tr, 


= lim 

n—too 


E L ^ w 

\v\=n 


lim Z n w a.s. and in C 8 . 

n—too 


Hence, Z = Z w and thus w ^ 0. Further, taking expectations in (3.39), we 
obtain 


= E [Z] 


E 


Y T Az] 


3 


L i>! 


w 


K[Zi\w. 
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i.e., w is an eigenvector corresponding to the eigenvalue 1 of E[Zi]. 

Conversely, suppose that w = E [Z] exists and is an eigenvector correspond¬ 
ing to the eigenvalue 1 of E[Zi] and that (Z n w) ne n 0 is uniformly integrable. 
Then Z n w —> Z w a.s. as n — > oo by the martingale convergence theorem and 
E[Z*"] = w. Since w is an eigenvector, we have w ^ 0 and thus P(Z W ^ 0) > 0. 
One can then check that Z w satisfies (3.39). 

To finish the proof in the case 1 < a < 2, we have to show that E[|Zi«;|^] < 
oo and m(f3) < 1 is sufficient for the /^-boundedness of (Z n w) ne ^ 0 . This, 
however, is a standard application of the Topchii-Vatutin inequality for mar¬ 
tingales. For details, we refer to [3, p. 182] and [71]. 

(d) Next, let a > 2. If w £ M. d is an eigenvector corresponding to the 
eigenvalue 1 of Z\ a.s., then Z n w = w a.s. for all n £ No and Z w = Z solves 
(3.39). For the converse implication, assume that Z solves (3.39) and that 
P (Z 7 ^ 0) > 0. By (3.41), Z £ C s for all s £ [l,a). Pick some s £ (1,2) if 
a = 2 and s = 2 if a > 2 and use the lower bound in the Burkholder-Davis- 
Gundy inequality [31, Theorem 11.3.1] to derive lower bounds for E[| Z — s ] 
as in the proof of Theorem 2.3 in [43]. It follows that E[| Z — tn| s ] = oo unless 
Z\W = w a.s. We refrain from providing more details. □ 

The proof for the case a. = 1 is much more involved, some parts of it will 
be deferred to the Appendix. As a preparation, we show that it constitutes no 
loss of generality to assume the stronger assumption (A5) instead of (A4). 

Lemma 3.6 Suppose that (A1)-(A3) and either (A4) or (A4’) hold. Let 4> 
be a solution to (3.2). Then there is a weight sequence T' := (Tj)jL 1 with 
N' < oo a.s. such that (A1)-(A4) resp. (A1)-(A4’) persist to hold for T', 
is a solution to equation (3.2) associated with T' as well and the martingale 
limit W defined in (1.8) is the same for T and T'. Moreover, ||Tj|| < 1 a.s. 
for all 1 < j < N' and the group G generated by T' is the same as the one 
generated by T. 7/T satisfies (A4’), then also the group O remains the same. 
If T satisfies (A4), then (A5) holds for T'. 

Proof We first deal with the case that (A4’) holds. Set T' := (T(v)) ve c( o) 
(recall that ||T(v)|| < 1 a.s. for all v € C(0)). Then [44, Proposition 3.7] states 
that if (T(n))|„| = i satisfies the assumptions (A1)-(A4’), then so does the family 
(T(v)) ve c(o)- That W is the same for the sequences T and T' follows from an 
application of [55, Theorem 9]. The additional point which we have to take 
care of here is that we need to know that the closed multiplicative subgroup 
of O (d) generated by the 0(v), v € C( 0 ), which we denote by 0 > , equals O. 
Clearly, 0 > C O. Conversely, let t > 0 and o £ O(rf) be such that to is in the 
support of Tj for some j € N. If we can show that o € 0 > , then it follows that 
O C 0 > . By (A2), there are t! £ (0,1) and o' £ O (d) such that t'o' is in the 
support of some Tk for some k £ N. Then o' £ 0 > by the definition of C(0) 
and 0 > . For the minimal m £ No such that t(t') m < 1, we have o(o') m £ 0 > , 
which implies o £ 0 > . 

Now we turn to the case that (A4) holds. Here, we obtain T' in two succes¬ 
sive steps. We start by considering as before the sequence T := (T(v)) ve c( o)- 
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As above, it follows that (referring to the proof of [44, Proposition 3.7] for 
the moment assumptions in (A4)) that (A1)-(A3) as well as the spread-out 
assumption and the moment assumptions in (A4) carry over. Also, as above, 
W is the same for T and T. That rj := E^,^^ ||r(t;)||“ £t(u)( - )] is spread- 
out can be obtained by arguments similar to those given in the proof of [ 6 , 
Lemma 1]. Hence there is an n £ N such that the n-fold convolution r]* n has a 
component which is continuous w.r.t. the Haar measure on § (d). Denote the 
corresponding density by / and assume without loss of generality that / has 
compact support and that c < f < d for constants 0 < c < d. Then if 2n has 
a component with density /* 2 , which is continuous [40, p. 295]. By the Fubini 
formula for the Haar measure on §(d) = R> x O (d) [33, Proposition 1.5.5], 

ri* 2n >jH 0{d)lD ®H R>]I (3.42) 

for some 7 > 0 and open sets D C O (d), I C R>. Replacing 2 n by An if 
neccessary, we may assume that D C SO (d). Since the latter group is compact 
and connected, we may invoke (the proof of) [16, Theorem 3] which gives that 
there is a k € N such that (H<o( d ^ D )* k > ei7so (d) for some e > 0. Then (3.42) 
yields that 

fj.(A,B) := rj* {2kn) (A x e~ B ), A C O (d),BC R measurable 
satisfies (M). On the other hand, 

^ = E 11^)11 <J (6( t ,)*-log(||T0)||))(') ’ 

■ \v\— 2 kn 

where (L(v)) v& y is the weighted branching process associated with T. It is 
readily checked that assumptions (A1)-(A3) and the moment assumptions of 
(A4) persist to hold for T' := {L(v))\ v \- 2 kn an d that also the martingale limit 
W is the same for T and T and thus the assertion follows. □ 

Remark 3.7 Under (A5), it holds that U = R> x SO (d) or U = R> x O (d). 
Since the group U' generated by T' will always be a subgroup of the one 
generated by T. it follows that under (A4), these are the only two possible 
cases for U. Note that U = R> x O(rf) while U' = R> x SO (d) is possible 
(for example, if ||Tj|| < 1 and det(Tj) = — 1 for all j = 1 This is 

not a problem since both the (U, a)- and (U', a)-stable laws are just the d- 
dimensional rotational invariant a-stable laws. 

Proof (Proof of Proposition 1.1 in the case a = A) By Lemma 3.6, we may 
assume throughout the proof, that ||T)|| < 1 a.s. for all 1 < j < N, and that 
(A5) holds if (A4) holds, i.e., we work with the sequence T' instead of T (but 
drop the superscript). 

We begin with the converse implication of (b) and suppose that w £ R d is 
an eigenvector corresponding to the eigenvalue 1 of E[Zi]. Then 

H = |E[ZrH < ||E[Z 1 ]|| • \w\ < E[ 11^11 ] • | w | < e[^ 11^11 j • | w | = \w\. 

L i> 1 J 
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Hence, equality must hold in this chain of inequalities and hence TjW = ||Tj|| w 
for all j £ N a.s. Consequently, Z n w = W n w Ww a.s. as n — > oo and Ww 
satishes (3.39). 

We are left with proving the direct implication in (b). Thus, assume that 
Z is ^-measurable and satisfies (3.39). First, we write M® = V+ + V+ where 
V+ = {x £ R d : ox = x for all o £ O}. Both, V + and V+ are invariant 

subspaces of for every o £ O. Write Z v + and Z v + for the orthogonal 
projection of Z onto V + and respectively. We have 

Z v ++Z v + =Z = ]T L(v)[Z\ v = Y L{v)[Z v +] v + Y L{v)[Z v +] v 

\v\—n |t?|=n |u| —n 

= Y \\L(v)\\[Z v +] v +Y L{v)[Z v t] v . 

\v\=n |i;|=n 

From linear independence (the left sum is in V+, the right sum is in Vq)~), we 
conclude that 

Z v + = 'Y / ||T(u)|| [Z v+ ] V a.s. for all n £ No 

\v\—n 

and Z y + = Y^ L{v)[Z v +\ v a.s. for all n £ No- 

\v\=n 

The equation for Z v + can be reduced to one-dimensional equations and it 
follows from [5, Theorem 4.13] that Z v + = Ww for some w £ V+. Taking 
expectations, we conclude that w = E[Z y +]. Now notice that V + = E\{¥\Z\\) 
and thus w is an eigenvector of E[2Ti] or zero. 

Observe that if SO (d) C O, w has to be zero: Using the change of measure 
(3.4) with fk{Li) := (0±w)k for 1 < k < d, we infer that 

W = E J2\\ T j\\(°j w ) =JE[OH> 

L i> i 

which implies w = 0\w a.s. Hence, SO(d) C O implies w = 0. Referring to 
Remark 3.7, we see that (A4) particularly implies SO (d) C O. 

It remains to show that Z v + = 0 a.s. We drop the V+ in the superscript 
and write Z for Z v + . We will use a variant of (3.40). To formulate it, consider 
the coming generation at time t > 0, which is C(t), dehned by Eq. (3.6). The 
idea is to switch from genealogical generations in the branching process to 
particles living roughly at the same time (with S(v) interpreted as the time 
of birth of particle v). The gain is that there is a better control over the birth 
times S(v) when v ranges over C{t) than when it ranges of {|f | = n}. Instead 
of considering the multiplicative martingales (M„(x))„ 6 n 0 fo 1 ' Z (defined via 
the characteristic function (j) of Z), we consider 

M C {t){x) = II cj)(L{v) J x), xeR d . 

v£C(t) 
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It can be checked along the lines of [2, Lemma 8.7(b)] and [5, Lemma 4.4] that 
Mc(t){x) —> M(x) for all x € l d with M(x ) = lim^oo M n {x) a.s. Arguing as 
in the proof of [43, Lemma 3.6] gives 

Z h = lim y E[L(v)[Z] v -\L(v)[Z] v \<h\T c{ t)\ a.s. (3.43) 

t—¥ OO z ' 

v€C(t) 

for every h > 0 with v({\x\ = h}) = 0 a.s. Here, Fc(t) is the cr-field that makes 
everything measurable in the weighted branching model defined by (C, T) that 
is born up to and including time t, see the proof of Lemma 8.7 in [2] for a 
rigorous definition. In the given situation, v = 0 a.s. and hence Z = Z h and 
h > 0 can be chosen arbitrarily, hence h = 1 is the most convenient choice. 
The key equation for us thus is 

Z= lim E L^/dlLHir 1 ) (3.44) 

£->-oo z ' 
v£C(t) 

= t h“ ( L ( v ) I ( et )+ i(«)( J T(||i(w)||- 1 ) - Ae*))) a.s. 

°° ■ vec(t ) vec(t) ' 

where I(t) := E[Z; |Z| < t], t > 0. For 0 < s < t, we have 

\I(t)-I(s)\ < E[|Z|;a < \Z\ < t] 

¥(\Z\ >x)dx- t¥(\Z\ >t) + s¥{\Z\ > s ) 

< f ¥(\Z\ > x) dx + s¥(\Z\ > s). 

J S 

Using this for the last summand in the last line of (3.44), we infer 



vec(t) 

^ /■luwir 1 _ 

= E / ¥{\Z\>x)dx+ E \\L(v)\\e t ¥(\Z\>e t ). 

v€:C{t) e v€C(t) 

(3.45) 


The second sum tends to 0 a.s. as t —> oo by (3.41) and the fact that 
Suec(t) l|£( w )ll = E[W|Jc( t )] —> W a.s. Regarding the first sum, using (3.41) 
and ||L(w)|| = exp(— S(v)), we obtain that, for arbitrary e > 0 , 



F(|^| > a:)dor < e E 11^)11 (S(v)-t) 

vGC(t) 


for all sufficiently large t. By the law of large numbers for (single-type) general 
branching processes [ 66 , Theorem 3.1], the last sum converges to a constant 
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multiple of W, see [5, p. 191] for a detailed argument. Since e > 0 was arbitrary, 
it remains to show that in (3.44) 

L(v)I(e t ) —> 0 in probability as t —» oo. (3.46) 

■uEC(t) 

Now notice that, for all t > 1, J(e 4 ) £ V+ and, by (3.41), 

r 5 * 

|J r (e*)j < 1 + J f‘(\Z\>x)dx < const • t. 

Hence, (3.46) follows from Lemma 3.9. □ 

Let us note one consequence out of the last step of the proof, which con¬ 
stitutes the analogue of [43, Lemma 4.9]: 

Lemma 3.8 Let X be a solution to (1.7) and let W h be defined as in (3.13). 
For each h > 0, there is a finite constant K > 0 such that 

\W h \ < KW a.s. (3.47) 

Proof Fix h > 0. It can be checked that the identity (3.43) holds with Z h 
replaced by W h . The subsequent estimates leading to Eq. (3.45) remain valid 
(also without restricting to V+). Subsequently, one has to use the tail estimate 
(3.24) rather than (3.25) and hence obtains that 


W h 


lim ^ L(^)/(e*) 

t—fOO Z ' 


< K'W a.s. 


for some K' > 0. If now lirnsup^^ |/(e 4 )| = oo, then \W h \ = oo on the set of 
survival, which yields a contradiction. Hence limsup^^ |/(e*)| < oo and the 
assertion follows. 


Lemma 3.9 Assume that (Al)-(A3), ||T)|| < 1 a.s. for all 1 < j < N and 
that either (A4’) or (A5) hold. Let V+ = {x £ : ox = x for all o £ 0} be 

the space of O -invariant vectors and let V = V+. Further, denote by (xt)t>o 
a bounded sequence in V. Then 


lim t 

t—foo 


Y, L fr) 


Xt 


0 in probability as t oo. 


(3.48) 


Proof Without loss of generality, we assume that sup t>0 |xt| < 1. 

Next notice that V+ is an O-invariant subspace of K d . Thus also its orthog¬ 
onal complement V is an O-invariant subspace of M. d . 

Step 1: First assume that (A4’) holds. In particular, O is finite. Let oi,... ,o p 
denote the elements of O with o\ = lo(d)- Then ((S(v), O(u))|„| =n ) n >o defines 
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a multi-type branching random walk with type space O. Pick an arbitrary 
x £ V. Then 

p 

= 0. (3.49) 

3 =1 

Indeed, for o € O, oO = O and, hence, o(oia; + ... + o p x ) = o\x + ... + o p x. 
Consequently, o±x + ... + o p x € V + . On the other hand, 0 \X + ... + o p x € V 
since V is O-invariant. This implies (3.49), which in turn implies 


W ^ 

- > ^ °3 X t 

p u 


0 


for every t > 0. Consequently, 


ueC(t) 


ty( y e- s ^o(v 

w 

- °3 

P 

i j= i 

- ' v£C(t):0(v)=Oj 

p 

y* 

( E e- 8 lv) -j)ot« 

v£C(t):0{v)=Oj 

3= 1 

P 

y* 

Y' e -s(v) _ y 

v£C(t):0(v)=Oj 


3=1 



The result now follows from [44, Theorem 2.14] once it has been checked that 
the assumptions of the cited theorem are satisfied in the present situation. 
The kind of checking that is needed to verify the assumptions can be found 
in the proof of [43, Lemma 3.14], it is here where the additional assumptions 
that 117} || < 1 a.s., 1 < j < N, enters. 

Step 2: Next, we turn to the more delicate case where (A5) holds, hence O is 
infinite. Suppose that I t .i, - - •, 7 tjPt is a partition of O and, for j = 1,. .. ,p t , 
define Oj := h^j f If o Ho(do) where Hq is the norrned Haar measure on O and 
htj = H<o(I t j). Due to the translation invariance of the Haar measure Hq, the 
matrix ht,\ 0 \ + ... + ht lPt o Pt is invariant under multiplication with elements 
from O and hence (h t pO\ + ... + h t>Pt o Pt )x = 0 for all x € V. Therefore, 


t 'y L{v)x t 

v£C(t) 


= t 


< t 


+ t 


y L(v)x t ~wy htjojXt 

vec(t ) j =i 

pt 

y y e~ s{v) ( 0 (v) -Oj)x t 

3 —1 v€C(t):0(v)El t j 

y ( e ~ s(v) -htjViAojXt 

3=1 v veC(t):0(v)eh,j ' 


(3.50) 
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Now assume that the partition is uniformly fine in the sense that 

max sup ||o- Oj|| = o(f _1 ) as t —> oo. (3.51) 

l =1 >-,pt oeh,j 

Since O C O(d) and O (d) is a d(d — l)/2-dimensional smooth manifold, it 
has box-counting dimension d{d — l)/2, see e.g. [36, Section 3.2]. Hence, for 
0 < e < | A1, with <5 as in (A4), we can choose the family of partitions in such 
a way that p t = o(t d< ' d ~ 1 ^ 2+e ) = o(t e ~ 1+e ) as t —> oo, this will be assumed as 
well. 

Then the summand in the second line of (3.50) vanishes as t —> oo. We 
thus centre our attention on the summand in the third line. We will prove that 

pt 

*£ £ e s( ' v ' 1 — htjWl —> 0 in probability as t oo. (3.52) 

!=i vec{t).o{v)eh,o 

We reformulate this problem in terms of a general branching process with type 
space O (see [45] for definition and details). For a measurable set ICO define 

<t>i(t) = e t ^e _Sj(0) l{o j (0) e /}l[o,Sj(0))( i )) t> 0. 

1> i 

For c > 0, let < j>j(t ) := l[o, c ] (!)</>/(!) and notice that 

Wi\v{t) = i[o,d (t)e* ^ e~ Sl ^ 1 fo( v )Oi (v)ef> Ifo.Sj(«)) (t), t > 0. 
l>i 

The general branching process counted with characteristic cj>f is defined as 
Z^{t):=Y J Wi\v(t-S{v)) = e t Y e- s Wl {OW£j} , t > 0. 

■hGV •yGC(t):5(v|| l; |_i)>t—c 

We write mf 1 := e”*E [Z^\, which is finite since it is bounded by E[.Z^°] = 1, 
see (3.7), and define 

h - : = (3 ’ 53) 
For 0 < c < t, we can then rewrite the left-hand side of (3.52) as follows: 


Pt 

tJ2 £ e-W-htjW 

Pt 

= tY {t) - h t , 3 W 

1 = 1 vGC(t):0(v)£l t ,j 

i=i 

p t . 

\ Pt 

< t Y (t) - e~ t Z 4,I *.i (t) 

)+tY( h ^~ hC iJ w ( 3 - 54 ) 

j= 1 

Pt 

+ tY e-'Z+kift-h^W . 

3= 1 

j=i 
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The first sum in (3.54) tends to 0 in C 1 as t —> oo when c = i/4. Indeed, 
observe that O = IJ^ implies that (jf It x (i) + ... + <f) c It (i) = t > 0 

and, hence, 


pt 


(3.55) 


^Z 0 J ‘o(i) = (t), t> 0 . 

i=i 

Then, taking expectations in the first sum in (3.54) and using Eq. (3.7) as well 
as E[W] = 1 leads to 


pt 


pt 


-t 


t E 

= t E 


e ^A»(f)-e-^A/t) 
i=i i=i 


E e SW (1 - l{ S - ( „||^|_ l) > t _ c }) 

€C(t) J 

= i P(SV( t )_i < t — c) < iP(5 T ( t ) — 5 r (f). i > c) —>■ 0 as t —> oo 

by Lemma C.l if c grows linearly with i. From now on, we fix c = i/4. The 
expectation of the second sum in (3.54) is bounded by 


E[Si] - E[Si A c] 
E[5r] 


E[5i 


/•OO 

t / P(Si > s) ds —> 0 

Jt/4 


as t —> oo since E[S'J] < oo (which implies P(S i > s) = o(s~ 2 )). 

The remaining sum in (3.54) requires much more attention. We begin the 
proof by observing that ].„(i — S(v )) = 0 for all v with S(v) < | due to 
the convention c = i/4. Consequently, for j = 1,... ,p t , 

e-‘Z^M(i) = ^ e- s ^ v) e-^- s{v)) [Z^] v (t-S{v)). 

vec (|) 

Using this, we can estimate as follows: 

Pt 1 

—t 


*E 

i=i 




pt 


„-S(v) 


,-S(v) 


m 


t—S(v) u Itj 




^ *E E e 

i =1 ^gC(|) 

pt 

+t E E e 

J—1 vGC(-|):5(f)<at 
Pt 

+*E E e_sw 

J=1 vGC(-|):5(f)>at 
Pt 4 

-'!>/„ E = : E J *w 

■? =1 wec(j) fc=i 
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for some (fixed) a € (|, 1). We consider each term separately. Using (3.53), the 
sum over (h c It .)j in Ji(t) is uniformly bounded by one, and [44, Proposition 
4.3] gives that linp^oo | Yl v eC( f) e ~ S ^ — VU| = 0 a.s. For J 3 (t), we change the 
order of the summation (the number of summands being finite a.s.) and use 
(3.55) and the change of measure (3.7) to obtain 


E[J 3 (t)] < tE Y e ~ SM J2( m t-sU+ h iJ 

vEC(^):S(v)>at j— 1 

< 2tP(S T(t/2 ) -i/2> (a-l/2)t), 

which vanishes by [44, Lennna A.3]. Turning to using that pt = o(t e ~ 1+e ) 
and EE„ e c(i) e -5 ^] = 1 (by (3.7)), we have that 


E[J 2 (<)] = tE 


pt 


-S(v) I m _ J,c 

/ , \ m t-S{v ) n It, 
3 =1 


E 

eC(±):S(v)<at 
Pt 

E l 

| m s t ' 3 — h c It . 

s>(a—\/2)t 

, ( s \ e + e | 4> c i ts , c 

< sup -- max \m s - h u 

s>(a-l/2)t — 1/2/ i<7<p t ' 


(3.56) 


To proceed further, we need a rate-of-convergence result from Markov renewal 
theory. Indeed, applying the change of measure, 


<t>r 


= e -*E 


= e“‘E 


= Z E E' 

n—0 L |i>|=r 

with 


Y^tM 1 ~ 

5(u))^e“ s<(, ' ) l/ tij (0(OT))l [ o ) s i („)](t- S(v)) 

oo 

YftA°n>t-S n ) 


vev 


i> 1 


— S(v) rc 


- S{v)) 


= E 


n =0 


: = l[o ) c](r)E[l/ t , i (oOi)l[o,s 1 )(r)] < l [ 0 ,oo)(d p (‘ S 'i > r ) ='■ ffW- 

By (A4), the function g satisfies g(r) = o(r^ +5+1 ). From (3.53), we conclude 

h c hj = / I ft,j(o,r)H 0 (do)dr 

and hence we can apply Proposition C.2 to deduce that the last term in (3.56) 
tends to zero as t —> oo. 

We finally consider J\{t) and proceed as in [44, pp. 735-736]. For fixed f, 
define 

Z vJ := e-(‘- s ( w ))[ \Z*ki] v (t-S(v)) 
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and similarly Z v with I t j replaced by O, i.e., Z v = Ylj =1 Conditioned 
upon jFc(t/ 2 )> the Z v , v € C{t/2) are independent. Let Z' v ^ := Z v jl^ Zv < e sMy % 

4>j ’ 

m.' V j := E [Z' v j], and J[(t) as Ji(t), but with Z v j and r m t _g/ v \ replaced by 
Z' v j and m' V j, respectively. On the set {Z v < e s( - v ) for all v € C(t/ 2)}, 


Mt) = + E e S{V) i. m t-s\v) - m v,j) 

i — 1 d£C(|) 

= + t e ~ S{V) ( m t-S(v)- m 'v)- 

vec (|) 

We want to prove that J\ ( t ) —> 0 in probability. To this end, we use the above 
decomposition and obtain for arbitrary r] > 0 , 

P(| > r/) = EfPdJ^t)! > 4^(1))] 

< El E nZv>e s(v) |^c(|))] + E[P(|J((t)| > vm^cd))] 
~ve e(f) 

2t~ 

—E 


V 


E e s(,,) ( m t-s(t,)- m 0 

«ec(|) 


The first and the last term can be dealt with as the corresponding terms in 
[44, pp. 735-736]. It remains to consider the middle term. 


^[n\m\>r,/2\F c(i) )\ 

Pt r , 

<J2 EP \ t E e- s ^\Z l vJ -m l V!j \> V /(2p t ) 

J=1 L ' w€C(|) 




Of) 


A 4p 2 f 2 

£ E^ e 


< 


3 =1 

4P(i 2 


rj2 


E 


E e“ 2S( " ) Var[Z; i |J- c(|) ] 

^ec(|) 

Pt 

E e_2S(,,) EE[(z;,) 2 i^(i)] 

y€C(|) 


< ^V-E 

T 


2+2 


< 


4;pit 
rf 

4ptt 2 


rf 


■E 


E 


E e~ 2S{ Mzll { z v <eS M} \Fc^)] 

ec( |) 

Z 2 

w(Zu) wE,) 1{z ’'- eS( ” )}|J ' c(|) 


^ e’ 2S <”)E 
weed) 


E 


„-so)_£ 


SO) 


^„2 


«ec(|) 

4 g£/2^2£+2e 


supE 


supE 


77 2 h 2 e+s{e t / 2 ) s >o 


ft. 2 f+5(e s ’0)) J 5>o 

W(e” s ^ 0 (s)) 




0 as t —> oo. 
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Here we used the triangular inequality, the independence of Z' v ■ and t\, 
Chebyshev’s inequality and the facts that t >->• t 2 /h 2 e+s(t) and t i-»- t/ti 2 e+s(t) 
are increasing and decreasing, respectively, for large t, t(\ogt) 2e+2e /h 2 e+s(t) —> 
0 as t —> oo, and the finiteness of the supremum, which follows from (A4) and 
is proved as in [43, Lemma 3.14]. □ 


3.10 Computing T 

In this section, we finish the determination of the ^-measurable Levy triplet 
(W 7 , S, v) of solutions of the homogeneous equation. In order to so, we will 
make use of some results proved in Section 4 below, which are postponed since 
their proofs do not require probabilistic tools like branching processes that are 
used in this section. 

As a by-product, we prove Proposition 1.6. We remind the reader of the 
definition of the functions 77 “, rjif, rj 1 and the vector 7 1 in Eqs. (1.19)-(1.22). 
pf, V 2 are defined in terms of a (U, o)-invariant Levy measure v a , i.e., satis¬ 
fying (3.29); p 1 and 7 1 are defined in terms of a measure p on § d_1 . Note that 
(A4’) implies G = R, hence U = { t® : t > 0 } x Co, see Proposition 4.1. Here 
t® := e( lnt ^, and we choose Q such that ||f^|| = t. 

Proposition 3.10 Assume (Al)-(A3), let <f> be a solution to (3.2) and let 
<P = exp(!f - ) be the limit of the multiplicative martingales, given by Proposition 
3.3. 

(a) Let a € (0,1). There is a (U, a)-invariant Levy measure v a such that a.s., 

&(x) = — W\x\ a r]i(x) + iW\x\ a ri%(x) Vie M d . (3.57) 

(b) Let a = l. 

(bl) Assume (A4) in addition. There is a c > 0 such that a.s., 

V(x) = -Wc|s| VieK d . (3.58) 

(b2) Assume (A4’) in addition. There is a z £ R d with Tj\ z = z and 

a finite Cu-invariant. measure p on satisfying f ( x , s)p(ds) = 0 for 
all x £ such that a.s., 

T(x) = iW(z,x)+ W(r/ 1 (x)+i('y 1 ,x)) (3.59) 

(c) Let a € (1,2). There is a (HJ, a)-invariant Levy measure v a such that a.s., 

3/{x) = i(Z, x) — W\x\ a r}^{x) + iW\x\ a rj 2 (x) \/x£R d . (3.60) 

(d) Let a = 2. Then there is a positive semi-definite dx d matrix E satisfying 
o T Eo for all o £ O and a z £ R d with z = X^'>i Tj z a.s., such that a.s., 

t T y T 

SP(x) = i {z,x) — W —Vcc€R d . (3.61) 
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(e) Let a > 2. Then there is a z € R d with z = Tj z a - s -, such that a.s., 

W{x) = i (z,x) (3.62) 

Proof By Proposition 3.3, L> is a Levy-Khintchine exponent (see (3.11)) with 
an ^-measurable Levy triplet (W 'By Lennna 3.4, v = Wv a.s. for a 
deterministic (U, a)-invariant Levy measure D. By Lemma 3.5, X = WE a.s. 
for a deterministic covariance matrix S satisfying S = oSo T for all o £ O. 
Moreover, v = 0 if a > 2 and E = 0 unless a = 2. The burden of the proof is 
to determine the random shift W'. We now consider the cases separately. 

(a) Let 0 < a < 1. Using the evaluation of the Levy integral in Lemma 
B.l, we obtain from (3.11) that for all x £ R d \ {0}, 

L r (x) = i(W',x) + W( - \x\ a r]±(x) + i\x\ a ij%(x) + \{^ a ,x)) 

= \{W' + 7 a VU, x) — W\x\ a iii(x) + iW\x\ a ri 2 {x), (3.63) 

for functions 77“, defined as in (1.19) and (1.20), resp., in terms of v a := v. 
Using that r]f(u T x) = r]f(x), j = 1, 2, and (3.17), we infer that for all n £ N 

T(x) = E ^UL{v) T x) 

\v\=n 

= iJ2( L (v)[W'+r w }v,x)- E \\L(v)f[W} v \x\ a vUx) 

|i>|=n |u|—n 

+ i E H L ( w )ir [W]v\x\ a ri 2 (x) a.s. (3.64) 

\v\=n 


Combining (3.63) and (3.64) and linear independence of 1 and i, we obtain 
(W'+^W.x) + W\x\ a ri%(x) 

= E (L(v)[W' +'y a W\ v ,x) + E [W] v \\L(v)f !«(*)■ (3.65) 

|u|=n |u| =n 

Dividing by |ar| and letting |a?j —> 00, we obtain 

(w' + rw,y) = '£(L{v)[iY' + 'fW\ v ,v) (3.66) 

\v\=n 


for all y £ S d_1 , hence W’ + ^ a W is a solution to (1.10), i.e., is an endogenous 
fixed point. Since a < 1, Proposition 1.1 yields that W' + 7 a W = 0 a.s. Then 
(3.57) follows from (3.63). 

(c) For 1 < a < 2, we start from Eq. (3.65). Dividing by |ar|, but considering 
\x\ —> 0 this time, we obtain the identity (3.66). Hence, W' + 7 Q W is an 
endogenous fixed point. Proposition 1.1 thus implies (3.60). 

(bl) Assumption (A4) implies U = R> x O with O = SO (d) or O = ©(d), 
see Remark 3.7. Thus, the r.v. with Levy triplet (0,0, v) is rotation invariant 
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and 1-stable. By [72, Theorem 14.14], its characteristic exponent equals — c|x| 
for some c > 0 and we obtain that a.s. 

\P(x) = \(W',x)— Wc\x\. 

Using (3.17) and considering real and imaginary part separately, we obtain 
that 

w = L ( v WU 

\v\=n 

and Proposition 1.1 yields that W' = 0 a.s. due to the assumption (A4). 

(b2) Let p be the spherical component of the Levy measure P, given by 
Proposition 4.3. We start by proving that / (x, s)p(ds) = 0 for all x E Ei(Q J ). 
We consider two different cases: 

1. Suppose x € Ei(Q t ) D .Ei(Cu) =: V\ i.e., u T x = ||u|| x for all u € U. 
For this case, we adjust the proof given of [43, Theorem 4.10 (bl)] to the 
present situation. Note that on V', (1.7) reduces to a fixed-point equation of 
a smoothing transformation with nonnegative scalar weights. Let ei,...,efc 
denote an orthonormal basis of V' and write Wj := (ej, W' + 'yW) where 7 
is as in (B.3). Further, let Sj := (ej,s ) for j = 1 and s E R d . Thus, 

using (B.3) and the linear independence of i and 1, we obtain for all r E R>, 
j = a.s. 

rWj — W— [ rsj log(|rsj|) p(ds) 

7T J S'*- 1 

= [ rsj P {ds) 

— — 7T /®d-i 

v=n v—n ^ 

-^Z\\ L { v )\\[ w ]v- f rs j log(l rs il) p(ds)- (3.67) 

Assuming for a contradiction that f sd ~i Sjp(ds) ^ 0 for some j, we choose 
r > 0 such that L d -i rsj log |rSjj p(ds) vanishes. Hence, upon dividing by r, 

Wj = ^||L(u)||[W j ] v -^||L(u)||log(|| J L(u)||)[W] v - [ s jP (ds). 

— — 7T /<sd-i 

v—n v=n ^ 

(3.68) 

By Lemma 3.8, Eqs. (3.12) and (3.13) and Lemma 3.4, we deduce that there 
is some K > 0 such that | W 3 \ < KW a.s., and consequently also |[W 7 -]„| < 
K\W] V a.s. for all v E V. This together with (3.68) yields that a.s. 

^||L(u)||log(||L(u)||)[W],- / s jP (ds) <2 KW. 

7 r /sd -1 

v=n 

The assumption / §d _i Sj p(ds) ^ 0 implies that the left-hand side tends to 00 
a.s. since lim^oo sup H=n ||L(u)|| = 0 a.s. by (1.5) and X)|„|=„ ll-^MII \W}v = 
W a.s. Contradiction! Hence, f sd _ 1 (x,s) p(ds) = 0 for all x G V'. 
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2. Suppose x G £a(Cu)' L . Since p is Co-invariant (see Proposition 4.3), 
it follows that so := / sp(ds) is Cu-invariant as well, thus 0 = (x,Sq) = 
f(x,s)p(ds). 

Hence, by 1. and 2. together, 

J{x, s) p(ds) = 0 for all x € Ei(Q J ). (3.69) 


By Corollary 4.2 below, the infinitely divisible law with Levy triplet (0, 0, v) 
is operator-stable with exponent Q. For such laws, [59, Theorem 13] shows 
that validity of (3.69) is equivalent to the existence of 7 1 (given by the fornula 
(1.22), cf. [59, Proposition 12 and Eq. (19)]) such that 

&(x) ■■= i(y 1 ,x)+r) 1 (x) 

satisfies = tiP(x) for all t > 0. Moreover, since t® commutes with 

Cu (see Prop. 4.1), the Co-invariance of p and the definition of 7 1 imply that 
7 1 is Cu-invariant, and so is rj 1 . Hence, E[{t^) T o T x) = t\P(x) for all t > 0, 
o € Cu, i.e., & is the exponent of a strictly (U, l)-stable law. 

Thus we can write 


V(x) = + WE(x), (3.70) 

and this equals, using (3.17) and the strict (U, Instability of 

\P(x) = i( ^ L(v)[W' — ^ X W] V ,x) + E nx)\\m\\[w) v 

|i;|=n \v\=n 

= i{J2 L ( v )i w ' - x) + WE(x). (3.71) 

\v\=n 


Substracting (3.71) from (3.70), we infer that W' — 
fixed point, hence equals zW for some z with E[JT>i 
1.1. All in all, 


^ X W is an endogenous 
Tj\z = z by Proposition 


\P(x) = zW + 'P(x). 


(d) a = 2 implies v = 0 a.s. and hence 


«P(a) = i{W',x) 

The claimed properties of £ are proved in Lemma 3.5. Using again the linear 
independence of 1 and i and (3.17), we deduce, since x € \ {0} is arbitrary, 

that 

w = £ L(v)[W'] v , 

\v\=n 

hence W' satisfies (1.10). By Proposition 1.1, either W' = 0 or W' = w for a 
deterministic mi/0 which satisfies w = ]Cj>i Tjiv a.s. 

(e) In this case, v and E vanish, W' can be identified as in (d). □ 
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Proof (Proof of Proposition 1.6) The notion of (U, a)-stability implies, using 
uniqueness of the Levy triple ( 7 , E, v a ), that v a satisfies (3.29) and that o T Eo 
for all o £ O. For a ^ 1, we can argue as before, using the identity ' V{u T x) = 
||u||“ T(x) and letting ||u|| 00 resp. ||«|| —>■ 0 to prove that only r]f 2 or E 
remain. 

If a = 1 and U = {t® : t > 0} x C, then a (U, l)-stable law is in particu¬ 
lar operator-stable with exponent Q (see Section 4.2), and necessarily of the 
form ( 7 , 0 , 1 /“), where / f(x)v a (dx) = f R> J sd _ ± f(t Q s)^ p(ds)dt. Then [59, 
Proposition 12 and Theorem 13] give that (3.69) is equivalent to the existence 
of 7 G such that (7, 0, 1 /“) is strictly operator-stable. In addition, if (3.69) 
holds, then (y^O, v a ) is strictly operator-stable, and 7 1 inherits C-invariance 
from p. Thus, ( , y 1 ,0,i' a ) is strictly ([/, l)-stable and if ( 7 , 0 , 1 /“) is strictly 
(U, l)-stable as well, then also (7 1 — 7, 0,0) is (U, l)-stable, which implies that 
z := 7 1 — 7 satisfies u T z = ||u|| 2 for all u € U. □ 

Propositions 1.6 and 3.10 together show that all solutions of the homoge¬ 
neous equation are of the form Z + Yw. It remains to solve the inhomogeneous 
equation. 


3.11 Proof of Theorem 1.5: The converse inclusion 

Since we have determined all solutions to the homogeneous equation in Propo¬ 
sition 3.10 above, we can now finish the proof of our main result by proving 
the converse inclusion in Theorem 1.5 (the direct inclusion has already been 
proved in Section 3.3). 

Proof (Proof of Theorem 1.5: The converse inclusion) Let X be a solution to 
(1.1) with characteristic function (j>. Write T’n(x) := fl |«1 =n^(-^( w ) Ta; )> an< i 
notice that the multiplicative martingale associated with f> takes the form 
M n (x) = exp(i(a:, W*)) ■ <S> n (x) (see (3.8)). The assumption that W* —> W* in 
probability implies 

<& n (x) —> M{x)/ exp(i(lF*, x)) =: x) 

in probability. Arguing as in the proof of [4, Theorem 4.2], it follows that 
ip(x) := E[^(x)] satisfies the functional equation (3.2) of the homogeneous 
smoothing transform and that T>{x) equals the limit of the multiplicative mar¬ 
tingale associated with ip{x), hence <P(x) = exp(lP(x)) with T(x) given by 
Proposition 3.10. We conclude that 

<j)(x) = E[M(x)] = E[exp(i(W*,a:))^(x)] = E[exp(i(W*, x) + ^(a;))]. 


□ 
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4 Matrices and measures invariant under actions of similarity 
groups 

In this section, we study the property of ( U , a)-stability in detail, for arbi¬ 
trary closed subgroups U C §(d). We start by describing the general structure 
of such groups. This will allow us to relate {U, a)-stable laws to operator 
semi-stable laws, and to characterize Levy measures and covariance matrices, 
satisfying the invariance properties (3.29) and (3.37), respectively. 


4.1 Structure of U and polar coordinates 

Let U $7 O (d) be a closed subgroup of the similarity group S (d), a particular 
case of which is U, the closed subgroup generated by the Tj, j = 1,..., N. 
We write G for the image of U under the group homomorphism ||w|| and 
distinguish between the discrete case G = r z for some 0 < r < 1, and the 
continuous case G = R>. As before, t® := and the right-hand side 

denotes the matrix exponential. 

Proposition 4.1 Let Cu := U D O (d). Then there is a subgroup Au C U 
isomorphic to G such that 

U ~ Au k Cu, 

in particular, every u £ U has a unique representation u = ac with a € Ajj, c £ 
Cu- Moreover, in the 

— discrete case: Au = {A n : n £ Z} for some A £ U with ||A|| = r, 

— continuous case: Au = {t® : t £ R>} for a d x d-matrix Q = Q' + cld, 
where Q' is skew symmetric and c ^ 0. Q can be chosen in such a way that 
Au and Cu commute. 

Proof The structure of U is given by Proposition C.l in [24], commutativity 
in the continuous case is proved in Proposition D.13 in [24]. It is proved there 
that (in the continuous case) Q is an element of the Lie Algebra of U, i.e., e® £ 
U C S(d). Any orthogonal matrix is the exponential of a skew symmmetric 
matrix, hence a similarity matrix u with ||u|| ^ 1 is the product of an orthognal 

matrix times a scalar multiple of the identity matrix, which we represent by 
e (lnc)I d _ □ 

The multiplication • : (Ajj x Cu) 2 —> Au k Cu is defined via the conjugation 
action induced on Cu by elements of Au- 

(ai, ci) • (a 2 , c 2 ) = (aia 2 , a^" 1 cia 2 c 2 ), (ai, cr), (a 2 , c 2 ) £ A v x Cu- 

Notice that a^ 1 cia 2 c 2 £ Cu since Cu is a normal subgroup of U. Further note 
that when all elements of Au commute with all elements of Cu, multiplication 
on Au « Cu simplifies to 


(ai,ci) • (a 2 ,c 2 ) = (a 1 a 2 ,c 1 c 2 ), (ai,ci), (a 2 ,c 2 ) £ A v x Cu 
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and hence, in this case, Au x Cjj = Ajj x Cu, that is, U is isomorphic to the 
direct product Au x Cu- 

According to the two cases considered, we introduce generalized polar co¬ 
ordinates as follows. Let 


Su 


{x £ R d : \x\ = 1} = § d_1 
{a: £ : r < |x| < 1} 


if G = R>, 

if G = r z , 0 < r < 1. 


(4.1) 


Then any x £ \ {0} has a unique representation x = as with s £ Su and 

a £ Au- Notice that in general, a is not a scalar. For example, in the setting 
of cyclic Polya urns, Section 2.1.2, 

A u = {t c : f £R>}, C v = {( k : 0 < fc < 6}, S v = § 1 , 
where £ = cos(27t/ 6) + isin(27r/6) is a primitive 6th root of unity. 


4.2 Operator (semi)stable laws 

An infinitely divisible law on R d with characteristic exponent is called (A, c)- 
operator semistable if there is a d x d-matrix A, b £ ]R d and c £ (0,1) such 
that 

<F(A t x) = c ] P(x) + \{x, b) for all x £ R d , (4.2) 

see [39, Definition 1.3.6]. It is called operator stable with exponent Q if there 
is a matrix Q and a mapping s i—> b(s) £ such that 

x P((t®) 1 x) = W(x) + i(x,b(t)} for all x £ R d , t > 0, (4.3) 

see [39, Definition 1.3.11]. The law is called strictly operator (semi)stable if 
6 = 0 or 6(f) = 0, respectively. Recalling the definition of (U, a)-stability in 
(1.12), and using the structure of U given in Proposition 4.1, we obtain the 
following Corollary. 

Corollary 4.2 Let rj be a (strictly) (JJ, a)-stable law. 

1. If Au = {A n : n £ Z} with ||A|| < 1, then (4.2) holds with c = ||A|| Q , i.e., 
rj is (A, ||A||“) (strictly) operator semistable. 

2. If Au = {e s ® : s £ R}, then (4.3) holds upon rescaling Q such that 
Ijt^H = t, i.e., rj is (strictly) operator stable with exponent Q. 

Notice that Cu did not play a role in the above considerations, therefore 
(U, a)-stability is more restrictive than operator (semi(stability. 
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4.3 Levy measures invariant under similarity transformations 

Using the generalized polar coordinates introduced in Section 4.1 above, we 
can now describe the structure of Levy measures satisfying (3.29). Of course, 
here U can be any closed subgroup of S(d). We write H Av for the Haar measures 
on An — G, which is the counting measure in the discrete case, and the image 
of d t/t under the map t t® in the continuous case. 


Proposition 4.3 Let is be a Levy measure on R d \ {0}. Then the following 
assertions are equivalent: 

(i) is satisfies (3.29). 

(ii) There is a Cn-invariant finite measure p on S\j such that for all D-integrable 
f : R d \ {0} -> R, 

f f(x)P(dx) = f f f(ax)\\a\\~ a p(dx)H AfJ (da). (4.4) 
J J A u J Su 

Notice that since An is not compact, the Haar measure H AfJ on An is unique 
up to a positive scaling constant only. In the discrete case, we stipulate that 
H Av is the counting measure. In the continuous case, we stipulate that H Av is 
such that the pushforward measure of H Au under the map An R>, a i-»- ||a|| 
is dt/t. 


Proof We consider the discrete and continuous case separately. 

In the discrete case, (3.29) yields that v(A~ 1 -) = ||H|| Q U Recall the defi¬ 
nition of Sip from (4.1). Setting p := D(- fl Su), we first observe that by (3.29), 
is is Cu-invariant and hence so is p. Further, [39, Theorem 1.4.5] implies that 
for any Borel set B C R d \ {0}, 


OO r. 

v{B) = £ \\A\rp((A n B)nS v ) = / \\a\\~ a p^ 1 B) n S v ) H Av (da) 


lair 


' Su 


1b( ax) p(dx) H Av (da). 


This proves (i) => (ii) in the discrete case. 

For the converse implication, we invoke [39, Theorem 1.4.4] which gives 
that each measure satisfying (4.4) is a Levy measure with D(A~ n -) = ||^4|| ,lQ u, 
i.e., P(a~L) = ||a|| Q v for all a £ An- Let U 9 u = aoc with a o € An, c G Cn 
and let B be a Borel subset of R d \{0}. In general, An and On do not commute, 
but since (apc)^ 1 is an element of U with norm ||(aoc) -1 || = ||ao|| 1 , there is 
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a d G Ciu such that (aoc) 1 = a 0 V 




v{c 1 a 0 1 B) = 


l B) 


. Then 

= \\a 0 \\ a Hc'- 1 B) 



= llaoir [ [ 

J An J Sn 

= ||aoc|| a P(B), 


l B (c'ax) ||a|| _a p(dx) H Av (da) 
l B (ac' a x) ||a|r Q p(dx) H Av (da) 
1 B {ax) ||a|r“ p{dx) H Av (da) 


where the c' a = a~ 1 c'a £ Cg. In the next-to-last line, the Cu-invariance of p 
was used. Thus (ii) =>• (i) is proved in the discrete case. 

Turning to the implication (i) => (ii) in the continuous case, choose Q in 
such a way that || a = t. Then (3.29) for (t^) -1 becomes P((t^) _1 •) = tD(-). 
Define 

p{B) := v({t Q x : x £ B, t > 1}), BC S^ 1 . (4.5) 

By Proposition 4.1, t® and C \j commute for every t > 0. Thus, we infer from 
Eqs. (4.5) and (3.29) that 


p{cB) = v(c ■ {t®x : x £ 1 , t > 1}) = p(B) 

for all c£ Cj and all Borel sets B C S d_1 , i.e., p is Cu-invariant. Moreover, 
[39, Theorem 1.4.12] gives that 


HQ 


I f tc(t®x)t 2 dt p(dx) 

J Sn J 0 

I I lc(t Q x)\\t Q \\ a p{dx) 

•J K.^> J jSirr 


dt 

T 


' An J Sn 


tc(ax) ||a|| a p(dx) H Av (da). 


(4.6) 


Here we used the particular scaling of Q and the fact that dt/t is the push- 
forward measure of H Av under a i-»- ||a||. Thus the implication (i) => (ii) is 
proved. 

For the converse implication, we use that each u £ U is of the form u = t®c 
for some t > 0 and c £ Cy. [39, Theorem 1.4.11] gives that if D satisfies (4.4), 
then v is a Levy measure satisfying v( y {t®)~ 1 Q = tD(C). Validity of (3.29) for 
all u £ U then follows as in the discrete case. □ 


4.4 Matrices invariant under orthogonal transformations 

In this section, we analyze the structure of a positive semi-definite dxd matrix 
V satisfying 

oSo T = 


£ for all o £ O. 


(4.7) 
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We defined O as the closed subgroup generated by the i, but it can 

be any closed subgroup of O (d). 

The main result of this section is the following proposition. To formulate 
it, we recall two notions. We say that a subspace V of R d is O-invariant if 
oV = V for all o £ O, and O-indecomposable if it does not contain any 
nontrivial O-invariant subspace. 

Proposition 4.4 Let £ be a positive semi-definite symmetric matrix, satis¬ 
fying (4.7). Then there is a decomposition 

R d = V + ®V-®V 1 ® ...®Vi (4.8) 

into O-invariant orthogonal subspaces with the following properties: 

(i) Every o £ O is the identity mapping on V+ and minus the identity on V-. 
V+ and V- are the maximal O-invariant subspaces with these properties. 
Further, V+ and V- are £-invariant subspaces and the restrictions £\v± 
are positive semi-definite symmetric matrices. 

(ii) For each i = 1,... ,1, V) is O-indecomposable and £-invariant, and £\y t is 
a nonnegative scalar multiple of the identity mapping on Vi. 

One ingredient in the proof of Proposition 4.4 is the following variant of 
Schur’s lemma, see e.g. [57, Corollary XVIII.6.2]. 

Lemma 4.5 Let F be a family of real d x d matrices, and suppose that {0} 
and M d are the only subspaces of R d that are invariant for each matrix in F. If 
£ is a symmetric matrix that commutes with every matrix in F, then £ = cld 
for a constant c £ M. 

For the proof of Proposition 4.4, we need another lemma. 

Lemma 4.6 Let £ be a positive semi-definite symmetric d x d matrix and 
o £ O (d) such that £o = o£. Let 

R d = V+®V-®V l ®---®V k (4.9) 

be a decomposition of R d into orthogonal o-invariant subspaces where V± = 
E±i(o) and, for each i = 1 ,...,k, Vi is a 2-dimensional, o-indecomposable 
subspace on which o acts as a rotation by an angle n ^ ft; £ (0, 2tt). 

Then V+,V-,V\,... ,Vk are £-invariant as well. 

Proof We have o -1 = o T since o £ O (d) and multiplication of £o = o£ 
with o T from the left and right yields that o T and £ commute. Hence so does 
a := (o + o T ). As symmetric matrices, a and £ are diagonalizable. This implies 
that the eigenspaces of a are 17-invariant and vice versa. To see this, pick an 
eigenvector v of a corresponding to the eigenvalue A and notice that 

X(£v) = £av = a(£v), 

i.e., £ maps the eigenspace E\{a) into itself. 
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The eigenvalues of a can easily be computed: Let b £ 0(d) be such that 
bob T is in normal form, i.e., 


bob T 


fid' 0 \ 

— Id" 

Re i 


V 0 ReJ 


for 2 x 2-rotation matrices Rg i , i = 1,..., k. Then bo T b T = ( bob T ) T is in normal 
form, too, and b(o + o T )b T is a diagonal matrix, with diagonal entries +2, —2 
and 2cos(0,;), i = 1,..., k. The corresponding eigenspaces are V+, V_ and Vj, 
i = 1,..., k if all Oi are distinct. In this case, the proof is complete. 

Suppose there is an * £ {1,..., k} such that 6i = 9j for some j ^ i. Then 
Vi © Vj c E 2co s(0i )(a). Further, 


oSVi = SoVi = SVu 


i.e., SVi is o-invariant. By the reasoning in the beginning of the proof, EVi C 
-^2 cos (6i)( a )- Hence EVi = {0} C V or V, = Vj for some j with 9i = 0j. 

We show that EVi = Vj for i ^ j is impossible. Choosing a basis of 
£ , 2 cos(e i )( a )> H l ias t° be a symmetric matrix w.r.t. this basis. In particular, 
it has to permute Vi and Vj. Hence Vi © Vj is 17-invariant. With respect to a 
joint basis of Vj and Vj , 





for a 2 x 2-matrix A with non-vanishing determinant. Thus, det(E\y.^y .) is 
negative, which violates the Hurwitz criterion for positive semi-definiteness. 

□ 


Proof (Proof of Proposition 4-4) For each o £ O there is an individual decom¬ 
position of the form (4.9), we denote its components by V±(o) etc. 

Defining V + := f) og0 V + (o) and V_ := PloeO Fl(o), we obtain E- and CD- 
invariant orthogonal subspaces, on which each o acts as the identity or minus 
the identity, respectively. V := (V+ © V-) -1 is E- and O-invariant as well. 

We consider the set 

V := {Vi © • • • © Vi : {0} ^ Vj C V is a E- and O-invariant subspace}. 

We do not distinguish between decompositions that consist of the same sub¬ 
spaces, but in a different order. V is non-empty since it contains V. The set 
V possesses a partial order a decomposition is larger than another one if 
the former is a refinement of the latter. Any totally ordered subset of V has 
at most d elements since every strict refinement decreases the dimension of at 
least one subspace by at least 1. Now pick a totally ordered subset of V with 
maximal number of elements and within this subset pick the largest element, 
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V\ ® • • • ® Vi, say. Clearly, V\ ® ■ ■ • © Vj is maximal with respect to We claim 
that V\,Vi are O-indecomposable. If not, then there is a Vj which contains 
two proper orthogonal subspaces that are invariant under every o £ O, and 
hence also ^-invariant by Lennna 4.6. This contradicts the maximality of the 
decomposition V\ © • • • © V) with respect to 

Consequently, we can decompose £ according to a maximal element of V. 
In particular, it suffices to solve 

0 |y £\y Oty — £jy for all O £ O 

separately for V £ {V+,V-,Vi ,..., V)}. 

The restriction of £ to V+ or V- can be any positive senri-definite sym¬ 
metric matrix, for conjugation by o\y ± is the identity mapping for all o £ O. 
To identify £\v t , where V) is an O-invariant and indecomposable subspace, we 
use Sclrur’s lemma (Lemma 4.5). This yields that each £\y i is a scale multiple 
of the identity on V), and since £\y i is positive semi-definite, the scaling factor 
is nonnegative. □ 


A The Choquet-Deny lemma 

Given a probability measure // on the similarity group S (d), let U be the closed subgroup 
generated by the support of //—if // is the step distribution of the associated multiplicative 
random walk (L n ) n ^ g , see Section 3.4, then U = U. 

Lemma A.l Let if : U — > R be measurable and bounded. If 

[ V>( ug)p(du ) = 'ip(g) (A.l) 

Ju 

for all u €. U, then ifc is constant /a-a.e. 

This is a consequence of [38, Theorem 3]. For the reader’s convenience, we state that 
theorem and show how the lemma can be derived from it. 

In the following, let G be a locally compact, separable and unimodular group. A prob¬ 
ability measure /x on G is called aperiodic , if the closed subgroup generated by the support 
of /x equals G. Write [G, G] for the commutator subgroup, i.e., the group generated by the 
commutators [a, b] := (ba)~ 1 ab, a,b E G, and [G,G] for its closure. Let H C G be a normal 
subgroup of G. Then G acts on H by conjugation (inner automorphisms), i.e., 

g.h := g~ 1 hg , g£G,h£H. 

For AC H write 

A g := {g.a : g E G, a E A}. 

The action of G on H is said to be compact if for each compact A C H \ {1q}, 1^ the 
unit element of G, A G is relatively compact, i.e., has compact closure. Then [38, Theorem 
3] reads as follows: 

Theorem A.2 Let /x be an aperiodic probability measure on G. If [G,G\ is Abelian or 
compact and if the action of G on [ G , G] is compact, then the only bounded, measurable 
functions -ip satisfying 

f 4>(ug) fi(du) = if{g) for all g G G 


are the p-almost everywhere constant functions. 
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Following the proof of [24, Theorem A.l], we show how Theorem A.2 applies to the 
situation here, i.e., G = U is the closed subgroup generated by the support of f.i. Then p is 
aperiodic on U by the very definition of U. Referring to Proposition 4.1, there is a closed 
subgroup Ajj of U , which is isomorphic to a closed subgroup of the multiplicative group M>, 
and a normal compact subgroup Cjj = U DO (d), such that U / Cjj ~ Ajj. The groups Ajj 
and Cjj (as a compact group, see [33, Theorem 1.4.1]) are unimodular and hence, by the 
Fubini formula for the Haar measure on U = Ajj Cu , [33, Proposition 1.5.5], U is unimodular 
as well. 

Clearly, the commutator subgroup of U is a subgroup of ©(d), hence its closure is 
compact. Moreover, for any compact AC [U,U]\ {Id}, A u is again a subset of ©(d) since 

u.[a,b] = u~ 1 [a, b]u = o~ 1 [a, b\o, 

where u = ||u|| o with o E Cu C ©(d). Hence, A u is relatively compact as a subset of a 
compact set. 


B Evaluating the Levy integrals 

In this section, we compute 

I{ x ) = if l ' X ’ y) - 1 ~ i ( x >2/) 1 [o,i](ll'l)) p ( d 2 /), xeR d (B.l) 

for a deterministic (U, a)-invariant Levy measure v, i.e. satisfying (3.29). 

Lemma B.l Let v be a deterministic Levy measure satisfying (3.29) for some 1 ^ a E 
(0,2), and define I(x) via (B.l). Then, for 0/r E R d , 

I(x) = — \x\ a r]f (#) + i\x\ a r)2 (x) + i(x, 7 a >, (B.2) 

with functions 77 “, rff; defined in (1.19) and (1.20), respectively, rjf and rff, are bounded 
real functions satisfying pJ{uJ x) = r}^(x) for all u £ U, rr E \ {0}, j = 1, 2, and 77 “ is 
nonnegative. The vector 7 “ satisfies 07 “ = 7 “ for all o E C\ j. 

Proof Fix 0 / a; E R d and notice that I(x) is finite since v is a Levy measure. Further, 
according to Proposition 4.3, there is a Cu-invariant finite measure p on S\j such that (4.4) 
holds. 

We set 

a , = f-/yl[o,i](MMd 2 /) = -f AfJ f Sv a 3 -l[o,i] (M) ||a|r“p(da:) ^(da) if a < 1 
7 1 — Sa v fsy aa: ^(l,oo) (l a3; l) ll a l| —a P(d^) Ha v (da), if a G ( 1 , 2 ). 

The asserted Cfj-invariance follows from the Cu-invariance of p. 

Recalling the definitions 

Vi(x) = f (1 - cos{(x,y)))u a (dy) 

V2( x ) = f (sin«rr, 2 /» - l{c,>i}(^,3/»^“(dy), 

(B.2) holds, and it remains to prove boundedness and invariance properties of 77 “, i = 1,2. 
Let u E U, then, using (3.29), 

Vf(.u T x) = 1 f (l-cos((x,uy)))v a {dy) 

= f (l-cos((x,y)))v a (dy) = y?(x), 

and the invariance of 77 J is proved along the same lines. This implies in particular that the 
continuous functions 77 “ are determined by their respective values on the (relative) compact 
set 5u, hence the asserted boundedness follows. □ 
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For cx = 1, we can compute a meaningful expression for rj 1 (x ) only for x G 
Notice that this implies x G E\(tQ) for all t > 0. Hence, using formula (4.6) for t>, we obtain 

V 1 ^) = f ad l (e 1<(t ° )Ta: ’ s> - 1 - i((t Q ) T a;,s>l { |( t Q)T s |<i } )t _2 dtp(ds) 

= Ld-i f R (e 1<te,s> - 1 ~ i(«a:,.s)l{ t <i})t _ 2 dtp(ds) 

= -[ + i-(^>s>log|(a;>s>| N )p(ds)+i(7,a:> (B.3) 

J S d ~ 1 V 7T ) 

for a suitable 7 G M d , see [72, Theorem 14.10] for details. 


C Auxiliary results from (Markov) renewal theory 

Lemma C.l Let (S n )neN 0 be a random walk with i.i.d. increments, So = 0, E[Si] > 0 
and E[(S']^) 2 ] < 00 . Let r(t) := inf{n G No : S n > t}, t > 0. Then, for every 0 < a < 1, 

£P(5 r ( t )_i < at) 0 as t —>• 00 . (C.l) 

Proof E[(5^”) 2 ] < 00 implies limt_>.oo t 2 P(5i > t) = 0. Further, it is known from standard 
random walk theory that limt —>.00 t~ 1 E[r(t)] = E[S'i] —1 . Consequently, setting, for n G N, 
A n ■= (So < t,.. ., S n - 2 < t , 5 n -i < at}, we have 

tP(5 r(t )_i < at) = t E P(A n n {s„+i >t}) < t E P(A„)P(Si > (1 - a)t) 

n>1 n>1 

< t 2 P(Si > (1 - a)t) • 7 E p ( t ^) - n ) 

t n> 1 

= t 2 P(Si > (1 — a)t) ■ — -1 — —¥ 0 as t —> 00 . 

□ 

^4 rate-of-convergence result in Markov renewal theory 

Throughout this section, let ((O n , S n ))neN 0 be a random walk on O (d) x M with increment 
law fi , say. The components O n and S n may be dependent. We assume that g, satisfies the 
minorization condition (M) and that 

E[Si] > 0 and E[|Si|^ +1+ ' 5 ] < oo (C.2) 

for some i > 0. Let O be the closed subgroup of 0(d) generated by the support of 0\. 

Proposition C.2 Let i > 0 and g : R — y [0, 00 ) be a measurable function that is decreasing 
on [0, 00 ), with g(t) = 0 for all t < 0 and limj—^oo t^+ 1 + e g(t) = 0 for some 0 < e < (<5 A 1). 
Then 

T _°°. 1 1 r 00 r 

lim sup t l+t E "22 f(On,t- Sn) - —— / / f(o,r)H 0 (do)dr = 0. 

t ^°°l/l<S L^0 - 1 E t S l] Jo JO 

Here and below, sup|^| <9 means the supremum over all measurable functions / : O x R —>• R 
satisfying sup oe <Q) |/(o, x)\ < g{x) for all x G M. 

The main ingredient in the proof will be the use of regeneration techniques for gen¬ 
eral state space Markov chains as developed in [9,67]. We sum up what is needed in the 
subsequent lemma. 
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Lemma C.3 There is a measurable space (12, Q) together with a family of probability mea¬ 
sures (P 0 ,r)oGO reM. and sequences of random variables ((M n ,R n )) n > q and ( r n )n>i 'with 

lPo,r(((M n ,i? n )) n >o C *) — + S n ))n>0 C ') (C.3) 

for all o E ©, r E R. Further, the following properties hold: 

(i) There is a filtration (£/n)n€N such that ((M n ,R n )) n > q is Markov adapted to (Gn)ne 
and {r n )n> l is a sequence of predictable (Q n )neN-stopping times, i.e., {r n = k} E Gk— l- 
(nj There are probability measures v on O and 77 on R, 77 having a bounded Lebesgue density, 
such that for alln> 1 and o E O, underP Q , ((M Tn+fc , R TTi+fc -R rTl _i))o<fc< Tri+1 _ TTl _i 

7 S independent of (Mo, Ro, • • •, M Tri — 1, R Tn —1) and has law P^(((M fc ,R fc ))^ = o C •)• 
(mj For eac/i bounded measurable function f : O —»• R, 


Tl-l 

E /( M «) 


= E 


i/®»,[n] / f(o)H 0 (do). 
Jo 


(C.4) 


(iv) There are C, A > 0 suc/i that P 0 ,r('7"i > ti) < Ce /or all o G O, r G I and n E N. 


Here and below, we use the shorthand P^®^ = f R P 0 ,r t'(do) ? 7 (dr), the same nota¬ 
tion for expectations, and sometimes omit the initial value Ro if if is irrelevant. 

Now we turn to the proof of Proposition C.2. 


Proof (Proof of Proposition C.2) In order to prove this result, we will combine methods 
from [ 8 ] and [ 68 ]. Below, we describe the steps of the proof and defer the technicalities to 
several lemmata. By splitting / into its positive and negative part, it suffices to consider 
nonnegative functions which are bounded by g. 

Let ((M n , R n ))n> 0 and ( r n)n>0 be as i n Lemma C.3. Define Vo := 0 and 


Vn ■— YARrt + i-1 Rr fc -l) — ^r n +l — 1 Rri-1- 

k =1 

Then, under each P 0>r , ( Vn) n >i 1S a random walk with i.i.d. increments and increment 
law P l/( g ) 7 ? (R ri _i E •) which is absolutely continuous since the law 77 of Ro is absolutely 
continuous. Since E[Si] > 0 and 

1 = P(Sn/n —>• E[5i] as n —> 00 ) = P/ d ,o(Rn/n —>• E[5i] as n —>• 00 ), 

we deduce that E^jjVi] = E l/ 0 7 ? [ri]E[ 5 i]. 

Let / be nonnegative and set 


m := E„ 


Tl-l 


x; f(M k ,t~R k ) 


Then we can proceed as in [ 8 , Section 4] to obtain 


^ ^ / (On, t — S n ) — E/ d? 0 ^ ^ / ( M n , t — R n ) 


= E 


id, 0 


Tl-l 


/ (M n , t — Rn ) 


(C.5) 




= E/,,0 


OO r r n + 1 f 

^ P 'y ^ /(Mk ? t — (Rfc — Rr n — 1) — Rr n — 1) j GTn — 1 

L n=l L k=r n 
ri-1 


y ' f {Mn , t Rn ) +E/ d? o f(t ~ Rti —1 — Vn) 


=: RiM + F2(t). 
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We show in Lemma C.4 that P Jte E\{t) tends to zero, uniformly over all / with |/| < g. We 
rewrite 


r r 

^2 if) = / Eiv<g )rj V fit -S-Vn) 

Jr L„_n 


^oC-Rri-l G ds) 


and use (C.4) to infer that (recall that / > 0) 


E„®„[Vi] X^ (r) 


dr = 


1 


1 


ri — 1 


T [ f(M k ,r - R k )dr 
k=0 

r i-! r 

Y / /(A4,r)dr 
u-n Jr 


E[Si] 

Then the claimed convergence rate holds, if 


E^0r7[Vi] 

4-r f f f(o,r)H 0 (do)dr. 
oil Jo Jo 


[ sup t e+e 

^v(g)ri 

l/l<s 



^u®r}Wl\ JO 


f(r) dr 


Y f(t -S-V n ) 

L n =0 

tends to 0 as t —»• oo. This result will be established in Lemma C.5. 


1 


?i d , o(-Rti-i G ds) 


Lemmata needed in the proof of Proposition C.2 

Proof (Proof of Lemma C.3) Observe that ((oO n ,r + S n ))ne N 0 1S indeed a Markov chain 
on O x R and that the increments of Sn — S n — l are independent of the past. Since /z satisfies 
the minorization condition (M), we have that for all o E SO(d), B E SO(d), 

P(oOi G B) = P(Oi G o-'B) > 1 H 0 {(o~ 1 B) n SO(d)) > ^H so(d) (B). (C.6) 

If O = SO(d), this shows that (oO n )neN is a Doeblin chain on O (see [65, Section 16.2] for 
the definition). 

If O contains elements with determinant —1 as well, then readily O = O (d), for SO(d) C 
O, and the product of two matrices with negative determinant has a positive determinant. 
Moreover, this necessitates P(det(Oi) = —1) > 0. Then, for all o E 0\S0(d) and B C SO(d), 

P(o0 2 € B) > f l { det(o')=- 1 }P(°l e ( ooT'B) P(Oi G do') 

> P(det(0!) = -1) ^H sow (B). (C.7) 

Thus, (oOn)nGN is a Doeblin chain in the case O = O(d), too. Its unique invariant probability 
measure is given by the normalized Haar measure Hq on 0(d). 

Again by the minorization condition for fi, it follows that there is an absolutely contin¬ 
uous measure rj(dx) := %i(x)dx such that 

P(oOi GA,SiGB) > ^H so(d) (A)r,(B) (C.8) 

for all o E SO(d) and all measurable A, B. 

(C.6) and (C.7) yield that there is some q < 1 such that, for all o E O, 

P(oO fc g SO(d) for k = 1,2) < q. (C.9) 

It follows that (oO n ) ne N is (SO(d),7/2, is, l)-recurrent in the sense of [8], with is := LTso(d)- 
Then, ((M n , Rn))n£N 0 can be constructed along the same lines as in [8, Section 3]: Under 
Po 5 r, let ((M n , R n ))ne N have the same transitions as ( oO n ,r + S n )ne hut whenever O n 
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enters SO(d), an independent B(l, 7/2)-distributed coin is flipped. If 1 shows up, then 
(M n _|_i, Rn+i — Rn) is generated according to v <S> rj, this event we call a regeneration ; if 0 
shows up, then (M n +i,i? n +i — R n ) is generated according to (1 — ^)~ 1 (P((M n Oi, Si) G 
•) — ^u<S>r )). Thus, the total transition probabilities are still equal to that of (oO n ,r-\-S n )ne N- 
Let to = 0 and r n be the nth regeneration time, see [9,67] for details. This gives Assertions 
1 and 2, while Assertion 3 is proved in [9, Theorem 6.1]. The construction together with 
(C.9) show that at least every third step, there is a uniform positive chance for regeneration, 
this yields Assertion 4. □ 

Lemma C.4 Let g be as in Proposition C.2 and ((M n , R n )) n >o as i n Lemma C.3. Then 


lim 

t —>oo 


sup t e+e+1 
\f\<9 


E/,,0 


- Tl-1 

- 71 = 0 



In particular, limt_>.oo t^+ e + 1 ^(t) = 0. Moreover, 

lim t* +e+1 Pr, 0 {R ri -l >t/2) = 0. 

t—* oo a ' 


0. 


(C.10) 


(C.ll) 


Proof In order to prove (C.10), we assume that g{ 0) < 1 and fix some |/| < g. Recall that, 
by Lemma C.3(iv), P/ d ,o( r i > n) < Ce~ Xn for some C, A > 0. Define nt = (log t)(i + 2)/A, 
t > 0. Then 




E /,,o 


J2 f(M„,t-Rn) 


< t l+t+1 E/,,0 


Tl-1 


T, 9(t - Rn) 


< t e+e+1 E /d ,o 


L n t J 

l{ri<n t } 9it - Ru)(1{r, 

71=0 


+ t: 


£+e+l 


E/,,o 


Tl-1 

l{Tl>71 t } 9{t~ Rn) 


<t/ 2} + l{H n >t/2})j 
= : Ii(t) + l2(t) Is(t). 


Recall that g is decreasing on [0, oo) and limt_>.oo t^ +1 + s g(t) = 0 for some 5 > 0. This gives 


r L™tJ 

h{t) = t e +* +1 E Id}0 l {T1 < nt} g(t-Rn)l{R n <t/ 2 } 

71=0 


< + He+1 (|mJ+l)g(t/2) -> 0. 


By Jensen’s inequality, E[|S n | K ] < n k 1 E[|5i| K ] for all k > 1. Thus, applying Markov’s 
inequality, 

l"*J L n tJ wr i<? K+i+«i 

nt) < t e+e+i x^ P (5„ >t/ 2) < ^ +e+i E 

o£+l+5 L n tJ £+1+5 

< ±- s —'£ f n i +*E[\S 1 \ i + 1 + 5 ] < 2^+ 1 + 5 ^-^E[|5i|^ +1+5 ], 

71=0 


which tends to zero as t -+ oo. Finally, 

h(t) < ^ +e+1 E, d , 0 [ril {T1> „ t} ] 


< t t+e+1 ^ntP/ d , 0 (n > nt) + J P/ d ,o(n > r) dr) 

< t t+e+1 (ii t Ce~ Xrlt + j e ~ Xnt ) 
t t+e+1 (n t + 1/A) 


< C 


tf+ 2 


as t —> oo. Thus we have proved the first and second assertion. (C.ll) follows from (C.10) 
with g(s) = l[ 0 ,t/ 2 )(s). □ 
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Lemma C.5 It holds that 


f sup t l+t 


l/l<s 

L 


X /(* -S~Vn) 


fp°° fir) dr 




P/ o(i? T1 _i e ds) = 0. 


(C.12) 


Proof We start by proving that the functions / are uniformly directly Riemann-integrable 
over |/| < g. Write 


C r i_1 

fit) = E u(8)6o £ f(M k ,t-s-R k ) 

J L fc=0 


77(ds), 


then observe that 


/ 




- ri -1 

Y f(M k ,t — R k ) 


L k =o 


d t < 


r T1_1 r 

X 9(t- Rk) 

L fc=o J 


di KqSq 



g(t)dt 


and recall that rj has a bounded Lebesgue density. Consequently, / is the convolution of a 
Lebesgue integrable function with a bounded function and hence continuous, see [7, Lemma 
VII. 1.2]. Further, arguing as in [ 8 , Eq. (4.4)], 

XI SU P l/(*)l < Ei/®jj[n] f X SU P l/(M)| #o(do) 

[lh,(l-\-l)h] J O tG [l2h,(l-\-l)2h] 

< Ew®jj[n] X SU P »(*) = : c 9- 

ZGZ [I2h, (Z+l)2/i] 

Hence, it suffices to show that g is directly Riemann-integrable. The latter is clear since g 
is monotone and Lebesgue-integrable (see [7, Proposition V.4.1(v)]). We have the uniform 
bound (cf. [7, Theorem V.2.4(iii)]) 


sup sup E u ®r, 
teR\f\<g 


X i/(* - y »)i 


L n =0 


X*i —4h,4h\ (Vn) 


n =0 


= : Dg < oo. 


Now we decompose the integral inside the limit in (C.12) according to the set {R Tl — \ > t/ 2} 
to obtain the following upper bound 


sup sup t^ +e 
s<t/2\f\<g 


E 




X /(* - s ~ V r) 

n =0 


f(r)dr 


HWVi] Jo 

+2D g t l + t W Idfi (R T1 _ 1 > i/2)- 


The second term tends to zero by (C.ll). For the first term, we invoke [ 68 , Theorem 4.2(ii)] 
(with G = <5 0 ), which gives (note that |/| < g implies |/| < g) 


lim t i+e 

t—> OO 


sup 

|/|<9 


OO 


V n ) 

n =0 


1 

EiWVi] 



0 


as soon as V\ has positive drift (here, E I/ 0 r? [Vi] = E^ 0 7 ? [ri]E[S'i] > 0 by the proof of 
Proposition C.2), a spread-out law (here, the law of Vi is even absolutely continuous) and 
^v<8>r) [|Vl |^+ e + 1 ] < oo (which is true by (C.2) and Lemma C.3(iv)) and g is bounded, 
Lebesgue-integrable and satisfies 

r2t 

t^ +e / g(r)dr —> 0 and t^ +e sup g(r) —»■ 0 as t —► oo. (C.13) 

J t r~>t 


Lemma C.4 gives limt^oo t^ eJtl g(t) = 0, which is sufficient for (C.13) to hold. 


□ 
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